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ABSTRACT 


In  recent  years,  a  number  of  advances  have  been 
made  in  determining  the  dynamic  properties  of  a  struc¬ 
ture.  These  advances  have  opened  the  field  of  modal 
analysis,  the  term  used  to  describe  the  use  of  experi¬ 
mental  frequency  response  data  to  determine  the  modal 
properties  of  natural  frequencies,  damping  ratios,  and 
mode  shapes.  One  of  the  main  problems  in  the  area  of 
modal  analysis  is  the  lack  of  complete  documentation 
for  the  techniques  in  use.  The  purpose  of  this  thesis 
is  to  investigate  some  of  the  modal  analysis  techniques 
in  use,  and  to  document,  fully,  the  assumptions  and 
the  equations  and  computer  algorithms  associated  with 
each  method.  Three  curve  fitting  techniques  are  pre¬ 
sented  along  with  some  example  problems  which  demon¬ 
strate  the  limitations  of  each  method 
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CHAPTER  I 


INTRODUCTION 


The  use  of  experimental  methods  to  determine 
the  dynamic  characteristics  of  a  system  or  structure 
has  become  widely  used  and  accepted  during  the  1970's, 
primarily  because  of  the  advances  in  digital  test 
equipment.  The  use  of  these  methods  has  become  com¬ 
monly  known  as  modal  analysis.  Experimental  modal 
analysis  provides  the  ability  to  pinpoint  vibration 
problems  by  showing  troublesome  frequencies,  damping 
factors,  and  mode  shape  data  which  can  be  used  by  a 
designer  to  pinpoint  weak  spots  in  a  structure. 

In  order  to  obtain  meaningful  answers  to  vi¬ 
bration  problems,  it  is  essential  that  an  adequate 
mathematical  model  be  selected.  If  a  model  is  avail¬ 
able,  there  are  methods  which  can  be  used  for  com¬ 
puting  the  dynamic  response  of  the  structure.  However, 
if  the  parameters,  i.e.,  mass,  stiffness,  and  damping, 
are  not  known,  or  if  the  model  itself  is  not  known, 
then  the  solutions  cannot  be  obtained.  One  approach 
to  modeling  is  through  the  use  of  transfer  functions. 

If  the  applicable  transfer  functions  can  be  determined. 
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then  the  response  to  any  input  can  be  determined. 

One  advantage  in  using  the  transfer  function  approach 
is  that  only  those  points  of  a  structure  which  are  of 
interest  need  be  considered.  The  determination  of 
these  transfer  functions  is  the  main  thrust  in  exper¬ 
imental  modal  analysis. 

In  general,  modal  analysis  involves  the  de¬ 
scription  of  the  dynamic  characteristics  of  a  struc¬ 
ture  by  observing  its  natural  modes  of  vibration  and 
the  associated  modal  properties,  i.e.,  frequency  and 
damping  ratio.  Natural  frequency  and  damping  ratio 
are  global  properties  of  the  system  and  consequently 
can  be  measured  at  almost  any  point  on  the  system. 

The  key  to  reliable  modal  parameter  estimates  is  the 
accurate  measurement  of  transfer  functions  of  the  sys¬ 
tem.  The  transfer  function  is  a  complex  valued  func¬ 
tion.  The  Laplace  transform  is  used  to  characterize 
the  transfer  function  because  it  can  be  used  to  con¬ 
vert  a  measured  frequency  response  function  into  a 
partial  fraction  expansion  with  each  partial  fraction 
containing  information  on  each  resonance  of  the  sys¬ 
tem  in  the  frequency  range  of  analysis.  The  frequency 
response  function,  which  is  the  ratio  of  the  output  to 
the  input  (assumed  here  to  be  the  ratio  of  output  dis- 


placement  to  input  force)  evaluated  along  the  jw  fre¬ 
quency  axis,  can  be  used  to  determine  the  frequency, 
damping,  magnitude,  and  phase  of  each  resonance. 

The  purpose  of  this  thesis  is  not  to  describe 
the  techniques  used  in  obtaining  the  frequency  re¬ 
sponse  function,  but  rather  to  explore  three  methods 
of  extracting  the  modal  parameters  from  the  measured 
data  and  to  point  out  the  advantages  and  disadvan¬ 
tages  of  each  method.  For  a  description  on  frequency 
response  measurement  techniques  in  mechanical  systems, 
the  reader  is  directed  to  references  (10)  and  (11).^" 

For  a  brief  review  of  the  literature,  prior 
to  1968,  concerning  the  modeling,  or  measuring,  of 
the  transfer  function,  the  reader  is  referred  to  a 
PhD  dissertation  by  William  R.  Shapton,  (15).  Shapton 

discusses  briefly  the  problems  with  earlier  attempts 

* 

at  transfer  function  measurement.  Some  of  the  earlier 
works  reviewed  by  Shapton  include  Savant  (14) , 

Chestnut  and  Mayer  (2),  Crafton  (3),  Tse,  Morse,  and 
Hinkle  (17) ,  Jacobsen  and  Ayre  (5) ,  Brooks  (1) , 

Simmons  (16),  and  many  others.  Shapton' s  work  was 

^From  this  point  on,  the  nomenclature  (  ) 

will  be  used  to  denote  a  reference  number  in  the 
bibliography. 


4 


in  the  development  of  procedures  for  predicting  the 
response  versus  time  of  surface  strain  of  a  complex 
structure  resulting  from  a  dynamic  excitation.  His 
assumptions  were  that  the  structure  was  capable  of 
being  represented  as  a  linear  system  with  time  in¬ 
variant  parameters,  and  that  the  system  could  be  ap¬ 
proximated  by  a  finite  number  of  ordinary  differential 
equations  with  constant  coefficients.  He  also  as- 
sinned  that  the  system  was  under damped  and  that  the 
damping  could  be  represented  as  an  equivalent  vis¬ 
cous  damping.  Shapton  presents  a  method  to  determine 
the  transfer  function  by  extracting  it  in  analytical 
form  from  the  experimentally  determined  frequency 
response  data. 

Most  of  the  works  performed  in  the  area  of 
deriving  transfer  functions  from  experimental  data 
have  as  their  initial  basis,  a  method  presented  by 
Kennedy  and  Pancu  (6) .  Kennedy  and  Pancu  presented 
one  of  the  first  papers  on  determining  modal  charac¬ 
teristics  from  test  data.  The  object  of  their  paper 
was  to  bring  out  an  awareness  of  the  phase  relation¬ 
ship  in  trying  to  correlate  between  measured  and 
calculated  modes.  They  presented  a  technique  to 
identify  normal  modes  from  polar  plots  (plot  of  re- 
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sponse  vector  with  forcing  frequency)  of  the  fre¬ 
quency  response  function.  The  basis  for  this  tech¬ 
nique  is  that  near  a  resonant  frequency  the  polar  plot 
will  approximately  describe  a  circle  (for  a  single 
mode) .  This  method  will  be  discussed  further  in  a 
later  section  of  this  thesis. 

A  more  recent  work  in  this  area  was  done  by 
Klosterman  (8) .  He  presented  new  methods  to  determine 
a  complex  modal  representation  from  experimental  data. 
He  generalized  the  method  of  Kennedy  and  Pancu  to 
include  the  case  where  damping  is  not  proportional  to 
the  stiffness  or  mass. 

Using  Klosterman' s  work,  more  recent  papers 
by  Richardson  (11)  and  Potter  and  Richardson  (9)  have 
presented  the  basis  for  measuring  transfer  functions 
using  a  least  squares  estimator  to  measured  frequency 
response  data. 

This  thesis  will  discuss  the  circle  fit 
method  of  Kennedy  and  Pancu,  and  two  other  curve  fit 
techniques,  one  which  uses  the  total  response  curve, 
and  one  which  uses  only  the  imaginary  part  of  the 
response.  Chapter  II  presents  the  dynamic  model 
used  for  an  elastic  structure  and  the  compliance 
transfer  function  for  the  structure  in  the  Laplace 


domain.  Chapter  III  discusses  four  methods  of  ex¬ 
tracting  the  desired  modal  parameters  from  the  fre¬ 
quency  response  function.  Chapter  IV  will  present 
two  example  problems  using  three  of  the  methods  dis¬ 
cussed  in  Chapter  III.  And  finally,  Chapter  V  is  the 
summary  and  conclusions  of  the  methods  discussed. 


CHAPTER  II 


THE  STRUCTURAL  DYNAMIC  MODEL 


2. 1  Assumptions 

In  most  problems  in  modal  analysis,  the  model 
is  assumed  to  be  described  by  a  system  of  n  simulta¬ 
neous  second  order  linear  differential  equations.  In 
the  time  domain  these  equations  are  given  by 

M  x  +  C  A  +  K  x  *  f  (t)  (2.1) 

where 

x  (t)  -  n-dimensional  displacement  vector 
f  (t)  »  n-dimensional  force  vector 
M  »  n  x  n  mass  matrix 
C  =»  n  x  n  damping  matrix 
K  »  n  x  n  stiffness  matrix 

For  the  purpose  of  this  thesis,  it  is  assumed  that 
the  matrices  M,  C,  and  K  are  real  valued  and  sym¬ 
metric.  The  energy  associated  with  the  mass  and 
stiffness  matrices  is  stored  and  can  always  be  re¬ 
covered,  however,  the  energy  associated  with  the 
damping  is  dissipated  and  is  lost  from  the  system. 


7 


8 


The  damping  mechanism  in  this  formulation  is  assumed 
to  be  viscous,  i.e.,  damping  force  is  proportional  to 
velocity.  In  practice,  most  mechanical  structures  ex¬ 
hibit  rather  complicated  damping  mechanisms,  however, 
Richardson  and  Potter  (13)  have  shown  that  all  of  the 
damping  mechanisms  can  be  modeled  in  terms  of  energy 
dissipation  by  an  equivalent  viscous  damper. 

2.2  The  Transfer  Matrix 

The  Laplace  transform  of  equation  (2.1)  yields 

M  [s2X(s)  -  sx(o)  -  x(o)J  +  C  [sX(s)  -  x(o)] 

+  K  X(s)  -  F(s)  (2.2) 

Assuming  that  the  initial  conditions  are  zero  (or  al¬ 
ternately  combine  them  with  F(s)),  the  transformed 
equations  become 

J^Ks2  +  Cs  +  K j  X(s)  -  F(s)  (2.3) 

Defining  B(s)  as  the  system  matrix, 

B(s)  -  Ms2  +  Cs  +  K  (2.4) 

then  equation  (2.3)  can  be  written  as 

B(s)  X(s)  -  F(s)  (2.5) 


m 
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B(s)  is  often  referred  to  as  the  dynamic  flexibility 
of  the  structure. 

Now  define  H(s)  as  the  inverse  of  B(s)  (as¬ 
suming  that  the  inverse  exists).  Then  equation  (2.5) 
becomes 


X(s)  -  H(s)  F(s)  (2.6) 

H(s)  is  the  n  x  n  transfer  matrix,  each  element  of 
which  is  a  transfer  function,  i.e.,  h^  is  the  trans¬ 
fer  function  which  relates  the  response  at  point  i  due 
to  an  excitation  at  point  j .  Because  of  reciprocity 
h^.  *  hj i  for  1  and  j •  Taking  a  closer  look  at 
each  element  of  H(s)  reveals  that  each  is  a  rational 
fraction  in  s ,  i.e., 

x  2n-2  ,  .  2n-3  , 

u^S  »  <  •  »  * 

hij - aet  8(a) - 

+  b2n-2s  +  b2n-l 

- dec  B(e> -  (2.7) 

where 

det  B(s)  *  determinant  of  B(s) 

The  roots  of  det  B(s)  are  called  the  poles  of  H(s)  and 
are  the  values  of  s  for  which  det  B(s)  ■  0. 

If  we  assume  that  the  poles  of  H(s)  are  dis- 
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tinct,  i.e.,  each  of  multiplicity  one,  and  that  all 
modes  are  underdamped,  each  element  of  H(s)  can  be 
expanded  into  a  partial  fraction  form  to  give 


where 

sk  *  k*"*1  pole  of  H(s) 
s*  *  complex  conjugate  of 

A^  =  complex  residue  matrix  of  H(s)  evaluated 
at  s  -  sk 

•fa 

*  complex  conjugate  of 

Each  A^  can  be  found  by  multiplying  H(s)  by  (s  -  s^) 
and  then  evaluating  the  result  at  s  *  s^,  as  long  as 
the  poles  are  distinct. 

2.3  Modes  of  Vibration 

The  resonant  poles  are  the  locations  where 
s  -  s^.  These  poles  can  be  expressed  as 

sk  "  *°k  +  jwk 

where 


*  the  modal  damping  coefficient  of  the 
kth  mode 

-  the  damped  natural  frequency  of  the 
kCk  mode 

The  undamped  natural  frequency  is  given  by 

°k  *  V°k2  +  “k2 

and  the  damping  ratio  is  given  by 


Each  complex  conjugate  pair  of  poles  (sk  and 
s^)  corresponds  to  a  mode  of  vibration  in  the  struc¬ 
ture.  Figure  2.1  shows  the  poles  for  a  single  degree 
of  freedom  system  as  viewed  looking  down  on  the  s- 
plane . 

Modal  vectors,  or  mode  shapes,  are  defined  as 
solutions  to  the  homogeneous  equation 

B(sk)  ^  «  0 

where 

u.  -  n- dimensional,  complex  valued  modal 


vector 


B (s^)  -  system  matrix  evaluated  at  the  pole 
location  s  ■  s^ 


Figure  2.1.  Poles  of  a  Single  Degree 
of  Freedom  System 

Potter  and  Richardson  (9)  have  shown  that, 
when  the  modal  vectors  are  defined  as  above,  the  res¬ 
idue  matrices  (A^)  can  be  written  as 

ik  '  akSkH£ 


where 
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*  complex  valued  scaler 
T 

■  an  n  x  n  complex  valued,  symmetric 
matrix 

In  other  words,  the  residue  matrix,  A^,  can  be  deter¬ 
mined,  to  within  a  multiplicative  constant  (complex 
valued),  by  a  mode  shape  vector  u^,  which  is  the  so¬ 
lution  to  the  homogeneous  equation  (2.9). 

Using  equation  (2.10),  the  transfer  matrix 
H(s)  can  now  be  written  in  the  form 

H(s)  -  ^  (2.11) 

k*  1 

or  it  can  be  written  as  a  summation  of  n-conjugate 
pairs 


n  r  t  *  * 

H(s)  -  V  SS&  +  (2.1: 

tZ_i  Is  -  sk  8  -  sk  J 

k  »  1 

Equations  (2.11)  and  (2.12)  are  the  forms  of  the 
transfer  function  from  which  all  the  desired  modal 
parameters  (frequency,  damping  ratio,  mode  shape)  can 
be  extracted. 


Potter  and  Richardson  (9) ,  and  also  (11) , 
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have  shown  that  because  of  the  form  of  the  matrix 
(as  defined  by  equation  (2.10)),  only  one  row  or  col¬ 
umn  of  the  transfer  matrix  need  be  measured  and  ana¬ 


lyzed,  since  all  other  rows  and  columns  contain  re¬ 
dundant  information.  For  example,  consider  a  two 
degree  of  freedom  system  with  two  sets  of  complex  con¬ 
jugate  pairs  of  poles  and  modal  vectors. 


Expanding  (2.12)  for  this  case,  the  transfer  matrix 
can  be  written  as 
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p— 

f—  *  * 

*  *  — 

u12u12 

u12u22 

u12u12 

u12u22 

- * 

- 

s  -  S2 

S  -  S2 

* 

s  -  S2 

S  -  S2 

+ 

*  * 

*  * 

U22U12 

u22u22 

u22u12 

u22u22 

- 9F 

- * 

s  -  S2 

S  -  S2 

S  -  S2 

s  -  S2 

(2.13) 


Note  that  in  each  row  and  column,  the  numerators  con¬ 
tain  the  same  modal  vector  multiplied  by  a  component 
of  itself.  The  modal  vectors  and  pole  locations  can 
be  identified  from  any  row  or  column  of  the  transfer 
matrix  except  those  which  have  components  which  are 
zero  (node  points) . 

In  practice,  the  pole  locations  sk  and  one 
row  or  column  of  the  residue  matrix  are  actually 
identified  from  one  row  or  column  of  the  measured 
transfer  function  data.  Suppose  that  in  the  two  de¬ 
gree  of  freedom  case  discussed  above  the  first  col¬ 
umn  of  the  residue  matrix  for  the  mode  (A^)  is 
measured.  This  column  can  be  expressed  as 


El 


alullull 


alullu21j 


(2.14) 
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2 . 4  Complex  Mode  Shapes 

The  previous  development  of  a  structural  dy¬ 
namic  model  required  that  the  damping  matrix  be  sym¬ 
metric  and  real  valued.  If  we  impose  no  further  re¬ 
strictions  on  the  damping  assumption,  then  the  modal 
vectors  will,  in  general,  be  complex  valued.  When  the 
modal  vectors  are  real  valued,  all  points  on  the  struc¬ 
ture  reach  their  minimum  or  maximum  values  simulta¬ 
neously.  This  is  not  the  case  if  the  modal  vectors 
are  complex.  Recall  from  equation  (2.8)  that  the  com¬ 
ponent  of  the  transfer  matrix  for  a  single  mode  k  can 
be  written 

A* 

a  (S) - (2.18) 

s  -  sk  s  -  sk 


For  each  complex  conjugate  pair  of  poles  there  is  a 
corresponding  complex  conjugate  pair  of  modal  vectors. 
The  ijttl  element  of  H^s)  is  (after  dropping  the  ij 
subscripts) 


Vs) 


(2.19) 


where 


vk  *  akukiukj ’ 


the  complex  residue  of  mode  k 
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th 

from  the  ij  frequency  response. 

The  inverse  Laplace  transform  of  equation  (2.19)  is  the 
impulse  response  of  mode  k  for  point  i  in  the  structure 
with  input  at  point  j .  If  only  mode  k  could  be  excited 
by  a  unit  impulse,  the  time  domain  response  would  be 


where 


XfcCt)  -  vk  exp(skt) 

•k  k 

+  vk  exp(skt) 

=  2  exp(-akt)  |^Re(vk)  cos  wkt 

-  Im(vk)  sin  u)ktj 

-  2  exp(-ckt)  |vk|  cos(ukt  +  ak) 

2 | vk |  "  Pea^  amplitude  of  the  impulse  response 
(twice  the  magnitude  of  the  complex 
residue) 

Im(vk) 

=  arctan  y »  which  is  the  phase 

k 

angle  of  the  complex  residue 


Notice  that  when  ak  *  0  or  180°  the  residue  vector  is 
real  valued.  It  is  this  phase  angle  which  determines 
the  phase  of  the  complex  mode  vector.  If  there  were 
no  damping  in  the  system  the  mode  vectors  are  real 
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valued  and  all  points  on  the  structure  reach  their 
minimum  or  maximum  displacement  simultaneously,  i.e., 
is  0°  or  180°.  When  is  an  angle  other  than  0° 
or  180°  then  the  node  lines  (lines  of  zero  displace¬ 
ment)  will  not  be  stationary  as  in  the  undamped  case. 

It  is  this  phase  angle  that  is  neglected  in  some 
modal  parameter  estimation  techniques.  In  the  tech¬ 
niques  described  in  this  thesis,  this  phase  angle  is 
accounted  for  and  does  not  present  major  problems 

2.5  Deriving  Mass,  Stiffness,  and  Damping 
Matrices  From  Frequency  Response  Data 

In  (9) ,  Potter  and  Richardson  have  shown  that 
the  mass,  stiffness,  and  damping  matrices  (M,  K,  and 
C) ,  and  hence  the  system  matrix  B,  can  readily  be  re¬ 
constructed  from  the  measured  modal  vectors.  The  sum¬ 
mation  in  equation  (2.11)  can  be  written  as 

H(s)  =  9  A'1^1  (2.10) 

which  is  an  n  x  n  matrix  where  the  columns  of  Q_  com¬ 
prise  the  u^  modal  vectors 


and  A~1  is  a  diagonal  matrix  containing  all  of  the  s 
dependence . 


Recall  from  equation  (2.4)  that 


B(s)  -  Ms2  +  Cs  +  K 
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and  evaluate  B(s)  at  s  ■  0,  then 
K  -  B (0)  -  H(O)'1 

One  can  see  that  H(0)  can  be  obtained  by  setting 
s  *  0  in 

K  *  H(O)"1  -  (0  ACO)"1©1)'1  (2.21) 


Thus  the  stiffness  matrix  can  be  derived  from  the  mea¬ 
sured  modal  vectors  u^  and  the  identified  a^  and  s^ 
complex  scalers. 


Recall  also  that  since  H  equals  B 


-1 


H  B  =  I 


(2.22) 


Differentiating  (2.22)  with  respect  to  s  gives 


hb'  +  h'b  =  0 


(2.23) 


and  differentiating  a  second  time  gives 

HB"  +  2hV  +  h"b  «  0  (2.24) 

• 

where  the  prime  denotes  differentiation  with  respect 
to  s.  If  we  then  differentiate  equation  (2.4)  with 
respect  to  s  we  obtain 

B*  -  2Ms  +  C  (2.25) 
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which  when  evaluated  at  s  *  0  yields  the  damping 
matrix. 

C  -  b' (0)  (2.26) 

Rearranging  equation  (2.23)  gives 

*  • 

B  =  -H  H  B 

or 

B*  =  -B  h'b 

which  when  evaluated  at  s  ■  0  yields  the  damping 
matrix. 

C  -  -B(0)  h'(0)  B(0)  -  -K  h' (0)  K 
or 

C  -  -K  ^0  (A”1) ' (0)6TJ  K  (2.27) 

Differentiating  equation  (2.25)  with  respect  to  s  gives 

M  -  h  b"  (2.28) 

Rearranging  equation  (2.24)  yields 

b"  -  -2  H'Vb’  -  H'V'b  (2.29) 

-1  " 

Using  the  relation  B  «  H  and  evaluating  B  at  s  -  0 


gives 


M  -  -B(0)h' (O)b' (0)  -  %  B(0)h"(0)B(0) 

-  -K  h' (0)C  -  %  K  h"(0)K 
or 

m  ■  -k  [e (a-1) ' (0)eT]  c 

-%  K  [0(A"1)"(O)eTJ  K  (2. 

Thus  by  defining  the  matrices 


0 


(A-1)" (0) 


2a, 


S1  \ 


2a 


2n 


*2n 


the  mass,  stiffness  and  damping  matrices  for  the  mea¬ 
sured  system  can  be  obtained  from  the  poles  s^,  sca¬ 
lars  a^,  and  modal  vectors  which  are  determined 
from  the  experimental  frequency  response  functions. 


2.6  Summary 

In  summary  then  it  has  been  shown  that  if  the 
transfer  matrix  of  a  structure  can  be  measured  (in 
the  Laplace  domain)  then  all  of  the  modes  of  vibra¬ 
tion  of  the  structure  can  be  characterized,  and  a 
matrix  model  developed  for  analysis .  The  problem 
that  remains  is  to  measure  the  transfer  function. 
Since  the  frequency  response  of  a,  structure  is  the 
transfer  function  evaluated  along  the  ju>  axis,  it 
then  contains  the  necessary  information  required  to 


model  the  transfer  function  (or  matrix) .  In  practice 
then,  the  frequency  response  functions  are  measured 
and  then  curve  fitting  is  performed  on  the  data  to 
curve  fit  into  the  s -plane  transfer  matrix.  Chapter 
III  will  discuss  some  of  the  different  curve  fitting 


methods . 


CHAPTER  III 


MODAL  PARAMETER  ESTIMATION  TECHNIQUES 

3 . 1  Parameter  Identification 

When  a  linear  structure  is  excited  by  a  broad 
band  input  many  of  its  modes  of  vibration  are  excited 
simultaneously.  Therefore,  because  of  the  linearity, 
the  transfer  functions  are  really  the  sums  of  the  res¬ 
onance  curves  for  each  of  the  modes  of  vibration  (see 
figure  3.1). 


Figure  3.1.  Multi  Degree  of  Freedom 
Transfer  Function 

At  any  given  frequency,  the  transfer  function  repre 
sents  the  sum  of  the  motion  of  all  the  modes  of  vi- 
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bration  which  have  been  excited.  The  dashed  lines  in 
Figure  3.1  represent  a  number  of  single  degree  of 
freedom  modes  which  when  added  together  form  the  multi 
degree  of  freedom  transfer  function.  The  degree  of 
overlap  (modal  overlap)  is  governed  by  the  amount  of 
damping  of  the  modes  and  their  frequency  separation. 
This  modal  overlap  is  caused  by  the  contribution  of 
the  tails  of  adjacent  modal  peaks.  Figure  3.2  shows 
a  transfer  function  which  is  lightly  damped  and  has 
sufficient  frequency  separation  so  that  there  is 
little  modal  overlap. 


Figure  3.3  illustrates  modes  with  heavy  damping  and 
very  little  frequency  separation.  If  the  damping  is 
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light  and  the  peaks  of  the  response  are  sufficiently 


separated,  then  the  response  data  in  the  vicinity  of 
each  peak  can  be  treated  as  if  it  were  a  single  de¬ 
gree  of  freedom  system.  However,  if  the  modal  over¬ 
lap  is  heavy,  i.e.,  the  "tails"  of  adjacent  modes  con¬ 
tribute  significantly  to  the  magnitude  of  the  mode 
being  investigated,  then  a  single  degree  of  freedom 
approach  will  not  work.  For  these  cases,  the  param¬ 
eters  of  the  modes  must  be  estimated  simultaneously 
through  a  multi  degree  of  freedom  approach. 

Once  the  natural  frequencies ,  damping  ratios , 
and  mode  shapes  are  determined,  an  analytical  expres¬ 
sion  for  the  Laplace  domain  transfer  function  can  be 
determined.  Once  the  analytical  expression  is  ob- 
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tained,  it  can  be  used  to  predict  the  response  of  the 
system  to  a  particular  type  of  input. 

3. 2  Single  Degree  of  Freedom  Techniques 

If  the  frequency  response  data  exhibits  light 
damping  and  the  peaks  are  separated  sufficiently  so 
that  very  little  modal  overlap  is  present  then  the 
following  techniques  may  be  considered. 

3.2.1  Quadrature  Response  Technique 

Recall  that  the  transfer  function  of  a  linear 
system  for  a  single  mode  k  can  be  written  as 

* 

h.  (s) - —  +  - (3.D 

s  -  sk  3  '  sk 

where  the  ak  is  the  complex  residue  at  the  kth  pole. 

If  we  let  s  =  ju,  the  frequency  response  function  can 
be  written 

* 

alc  ak 

\  (j«)  =  - - -  +  - 1 F 

ju>  -  sk  ju)  -  sk 

* 

»  ak  +  *k 

jw  +  ok  -  jwk  ju  +  crk  +  juk 


(3.2) 
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The  magnitude  of  a  typical  frequency  response  function 
is  shown  below  in  figure  3.4. 


Figure  3.4.  Frequency  Response  Curve  For 
A  Single  Degree  Of  Freedom 


In  the  neighborhood  of  m  =  the  majority  of  the 
resonance  curve  can  be  approximated  by 


\(» 


ak 

ju)  -  sk 


(3.3) 


If,  for  convenience,  we  remove  a  factor  of  from  the 
residue,  then  equation  (3.3)  can  be  written 


rk  rk 

^k^^  ”  2j  (jw  -  s,  )  ”  2 ((a),  -  u))  +  jo,  ) 


(3.4) 


which  can  be  written  in  real  and  imaginary  parts 


Real 


V,  (  \  „i\ 

■B  V 

rlk(»k  -  (o)  +  r2kok 

*5l 

(a>k  -  w)  +  ak 
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Imag 


r2k<“k 


-  ">  -  rlk°k 


(co,  -  0>)  +  ot 


where 


rk  ”  -  rlk  +  Jr2k 

Sk  =  "0k  +  ^“k 


(3.5) 


When  a)  *  the  above 

expressions  become 

Real 

hk(jw)] 

r2k 

r  t 

rllf 

Imag 

hk(jo))J 

■ 

(3.6) 


In  cases  of  light  damping  (or  proportional  damping) 
the  residues,  rk,  are  nearly  (exactly  for  proportional 
damping)  real  valued,  i.e.,  r2^  »  0.  When  this  as¬ 
sumption  is  valid  then  the  natural  frequency  and  mode 
shape  data  can  be  extracted  from  the  quadrature,  or 
imaginary  part,  of  h^(juj).  The  natural  frequency  is 
the  frequency  at  which  the  peak  in  the  quadrature  re¬ 
sponse  occurs ,  or  the  frequency  where  the  coincident 
(real)  part  of  h^(jw)  is  zero.  The  mode  shape  con¬ 
tribution  is  simply  the  magnitude  of  the  quadrature 
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response  at  the  resonant  frequency.  The  damping  ratio 
can  be  determined  by  the  familiar  half -power  (band¬ 
width)  method,  i.e., 

"u2  " 

*  “  *  [  (3'7) 

where  and  ^  are  the  frequencies  at  which  the  re¬ 

sponse  is  reduced  to  — — .  times  the  peak  value,  and 

wn  is  the  natural  frequency.  As  stated  before,  this 
method  is  only  good  for  cases  where  the  damping  is 
very  light.  One  can  also  get  into  problems  with  this 
method  if  the  modes  are  closely  spaced  in  frequency 
and  there  is  significant  influence  on  each  mode  due 
to  adjacent  peaks.  Because  of  the  limitations  of 
this  method  it  will  not  be  discussed  further  in  this 
thesis . 

3.2.2  Method  of  Kennedy  and  Pancu 

In  1947,  Charles  Kennedy  and  C.  D.  P.  Pancu, 
published  a  paper  entitled,  "Use  of  Vectors  in  Vibra¬ 
tion  Measurement  and  Analysis"  (6) .  This  paper  set  the 
stage  for  the  present  day  techniques  used  in  modal  a- 
nalysis.  The  main  purpose  of  their  paper  was  to  iden¬ 
tify  the  normal  modes  of  vibration  of  mechanical  sys- 
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terns.  They  assumed  that  the  damping  in  the  system  was 
proportional  to  the  displacement  and  90°  out  of  phase 
with  it.  They  also  assumed  that  the  off-resonant  vi¬ 
bration  is  constant  in  phase  and  magnitude  as  the  sys¬ 
tem  passed  through  a  resonance.  The  authors  made  use 
of  a  vector  response  plot,  more  commonly  called  a 
polar  plot  or  Nyquist  plot,  which  plotted  both  real 
and  imaginary  components  of  the  frequency  response 
data  vs.  frequency.  An  example  of  a  Nyquist  plot  is 
illustrated  in  Figure  3.5.  Their  basis  for  using  this 
type  of  plot  was  that  near  a  resonant  frequency,  the 
polar  plot  would  describe  a  circular  arc.  The  plot 
for  a  single  degree  of  freedom  system  is  shown  in 
Figure  3.6.  In  terms  of  the  analytical  transfer  func¬ 
tion, 


hk(“)  (0).  -  W)  +  jo. 


(3.8) 


where 


I 

ak 

rk  |  = 

T~ 

^magnitude  of  the  complex  residue 


a,  =  phase  angle  of  the  residue 


If 


1  and  -  0  then  it  can  be  shown  that 


I 


Figure  3.5.  Polar  Plot  Of  A  Frequency 
Response  Function 
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Re 


[»k]2  +  ^[“k] 


1  v  2 


16a; 


(3.9) 


which  is  a  circle  in  the  Nyquist  plane  with  a  radius 
of  and  centered  at  -j (z^- )  •  The  complex  residue 


Figure  3.6.  Modal  Resonance  In 
Nyquist  Plane 

vector,  r^»  merely  expands  the  radius  of  the  circle 
and  rotates  it  in  the  Nyquist  plane  clockwise  by  an 
angle  away  from  the  negative  imaginary  axis  (see 


Figure  3.7).  Kennedy  and  Pancu  only  considered  the 


Figure  3.7.  Complex  Modal  Resonance 
In  Nyquist  Plane 


cases  where  the  residue  vector  was  real  valued,  i.e., 


It  was  shown  by  Kennedy  and  Pancu  that  the 
natural  frequencies )  of  a  system  could  be  determined 
from  a  Nyquist  Plot  by  observing  where  the  arc  length 
of  the  plot  was  a  maximum  for  a  given  change  in  fre- 
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quency  (see  Figure  3.8).  The  arc  length  will  increase 


Imag 


Figure  3.8.  Determination  of  Natural 

Frequency  From  Nyquist  Plot 

until  (u  passes  through  a  resonant  frequency,  then  the 
length  will  decrease  and  later  will  increase  again  as 
another  resonant  frequency  is  approached.  The  reso¬ 
nant  frequencies  can  be  determined  by  picking  off  the 
points  where  the  arc  length  reaches  a  local  maximum. 

Once  the  natural  frequencies  are  found,  use 


is  made  of  the  fact  that  near  a  resonance  the  arc  ap¬ 
proximately  describes  a  circle.  A  circle  is  then  fit 
to  the  data  near  each  resonant  frequency. 
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One  method  that  can  be  used  to  circle  fit  the 
data  is  the  method  of  least  squares.  The  general  equa¬ 
tion  of  a  circle  is 

x^  +  y^  +  ax  +  by  +  c  *  0  (3.10) 

where  a,  b,  and  c  are  constants  related  to  the  location 
of  the  center  of  the  circle  and  the  radius . 


(3.11) 


Squaring  this  error  and  summing  this  result  over  all 
the  descrete  frequency  points  in  the  area  of  the  de¬ 
sired  fit  gives 


I 


(x2  +  y2  +  ax  +  by  +  c)2 

(3.12 


To  minimize  the  error  squared,  the  partial  derivatives 
of  (3.12)  with  respect  to  a,  b,  and  c  should  equal 
zero.  Thus  we  have 


n 


X 


2  >  (x2  +  y2  +  ax  +  by  +  c)x  *  0 


k  -  1 
n 


2  >  (x2  +  y2  +  ax  +  by  +  c)y  *  0  (3.13 


k  -  1 
n 


X 


2  ^  (x2  +  y2  +  ax  +  by  +  c)  =0 


k  -  1 


In  matrix  form  we  get  equations  (3.14).  Using 
Gaussian  elimination,  the  constants  a,  b,  and  c  can 
be  found  which  then  will  give  the  location  of  the 
center  of  the  circle  and  its  radius .  Once  the  circle 
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is  determined  it  can  be  treated  as  a  single  degree  of 


freedom  system  where  the  diameter  of  the  circle  can 
be  used  as  the  magnitude  of  a  component  of  the  mode 
shape  and  the  phase  angle  can  be  computed  knowing  the 
location  of  the  center.  Consider  the  circle  shown  in 
Figure  3.9.  The  coordinates  of  the  center  of  the 
circle  are  now  known.  Also,  the  coordinates  of  the 
point  where  u  equals  the  natural  frequency  are  known. 
The  angle  can  then  be  determined  from  the  following 
relation 


arc tan 


I*  -  *c 

|y  -  yc| 


(3.15) 


where 


x„  ■  Real  axis  coordinate  of  the  circle  center 
c 

yc  ■  Imaginary  axis  coordinate  of  the  circle 
center 

The  sign  of  the  angle  can  be  determined  from  the  direc 
tion  of  rotation  from  the  imaginary  axis .  A  positive 
angle  is  rotated  clockwise  from  the  ju>  axis. 


Figure  3.9.  Circle  Coordinated 


The  damping  ratio  can  be  determined  in  dif¬ 
ferent  ways .  One  way  provides  for  the  ratio  to  be 
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calculated  at  two  different  points  and  then  averaged. 
Consider  the  equation  of  motion  in  the  Laplace  domain 
of  a  single  degree  of  freedom  system 


(msz  +  cs  +  k)X(s)  =  F(s) 


The  transfer  function  can  be  written 


- 2 - T~ 

m(s  +  2^wns  +  (d^) 


<l/k)coJ 

- 1 

s4  +  2co)ns  +  wn 


or 


u> 


X(s)  =  ~n _ 

x<s>sta  ’  +  2!“n5  +  “n 


where 


2cu)n  =  c/m 


o)n  *  k/m 


X(s)  =  F(s)/k  =  "static  displacement" 

S  La 


Then 


H(ju>) 


X(jw)  _  <4 _ 

X(ju))sta  (w*  -  a)Z)  +  j2cu>wn 


(3.16) 


(3.17) 
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cj 


(oj^  -  hi1)2  +  (2ciou)n) 


j  e~i*  (3.18) 


where 


2?u)uj 

p  =  arc  tan  — ^ ^ 

u )  -  u) 

n 


(3.19) 


Looking  at  the  Nyquist  plot  for  this  function,  Figure 
3.10,  one  can  see  the  phase  angle  <j>  for  an  arbitrary  u. 


Imag 


Figure  3.10.  Nyquist  Plot  Of  Single  Degree 
Of  Freedom  System 


Now  define  0  =  2$  where  0  is  measured  from  the  diam¬ 
eter  of  the  circle  drawn  through  the  resonant  fre¬ 
quency.  Substituting  0=2$  into  equation  (3.19)  and 
solving  for  t,  gives 


(u>£  -  o)^)  tan  (0/2) 


2wu> 


(3.20 


n 


The  only  information  needed  to  get  the  angle  0  are  the 
coordinates  of  the  center  of  the  circle  and  the  coor¬ 
dinates  of  a  frequency  point  on  one  side  or  the  other 
of  the  resonant  frequency.  Consider  the  example  shown 
in  Figure  3.10. 


V  =  arc tan 


0  =  180 


X1  '  xc 


-  y. 


lyi 


Substituting  into  equation  (3.20)  gives 


/  2  2^  __  180  -  Y 

'wn  ’  ui'  tan  - 2 - 

2u, 

l  n 


The  same  thing  can  be  done  using  and  then  the  two 
values  of  c  averaged. 

This  method  of  circle  fitting  is  accomplished 
for  each  resonant  frequency  observed  in  the  Nyquist 


plot  for  a  given  frequency  response.  If  this  is  accom¬ 
plished  for  each  frequency  response  in  one  row  or  col¬ 
umn  of  the  transfer  matrix,  then  the  mode  shapes  can 
be  determined,  i.e. 


(al)n  *  diameter  of  circle  at  determined 
from  frequency  response  h^ 

(ai>2i  =  diameter  of  circle  at  determined 
from  frequency  response  h^ 


and  similarly  for  the  phase  angles  (a,)... 

K  1J 

Immediately  one  can  see  a  problem  with  this 
type  of  approach.  If  the  modes  happen  to  be  closely 
coupled  where  one  mode  dominates,  the  circle  fit  of 
the  smaller  mode  may  not  have  sufficient  points  to 
resolve  the  amplitude  of  the  residue,  i.e.,  the  di¬ 
ameter  of  the  circle  fit  through  the  smaller  mode  may 
be  incorrect  thus  giving  rise  to  errors  in  the  mode 
shape  and  also  the  damping  ratio  associated  with  that 
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In  order  to  reduce  the  limitations  imposed  by 
single  degree  of  freedom  assumptions,  one  must  consid¬ 
er  a  multi-degree  of  freedom  curve  fit,  which  is  the 
subject  of  Section  3.3. 

Chapter  IV  will  present  example  problems 
using  this  method  of  determining  modal  parameters. 
Also,  Appendix  A  gives  a  listing  of  a  computer  al¬ 
gorithm  derived  for  circle  fitting  frequency  response 
data. 

3 . 3  Multi  Degree  of  Freedom  Curve  Fitting 

For  cases  where  the  amount  of  modal  overlap 
is  sufficient  to  cause  significant  errors  by  use  of 
single  degree  of  freedom  techniques,  a  multi  degree 
of  freedom  technique  must  be  used.  Two  techniques 
will  be  described  in  this  section. 

3.3.1  Complex  Curve  Fit 

This ' technique  is  one  which  matches  the  sum¬ 
mation 


h(s) 


+ 


(3.21) 
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along  Che  frequency  axis  (s  =  ju)  to  the  measured 
frequency  response  function  data.  A  least  squares 
error  estimation  is  used  to  account  for  modes  that 
are  outside  of  the  frequency  range  of  interest. 

Consider  equation  (3.21)  evaluated  at  s  «  ju. 


n 

h(ju)  =  ^  - 


* 

*k 


k  =  1 


0)  -  s.  ju  -  s. 


(3.22) 


This  summation  is  over  n-modes .  The  summation  can 
be  rewritten  as 


h(ju>)  - 


m 


m 


j  U  “  s  ju>  -  sm 

m  m 


n 

v  A 

Z_l  JW  - 


k  =  1 
k  /  m 


ak 

_  ^ 
JU)  -  sk 


(3.23) 


If  we  define 


a 


m 


m 


ju  -  s 


m 


n 

y^ 

Z_J  ju»  - 


* 

*k 


k  -  1 
k  f  m 


s.  ju  -  Si 


(3.24) 


then  equation  (3.23)  can  be  written 


h(ju)  = 


jo,  -  sT 


am  +  j  (“  •  “W 


(3.25 


If  the  system  is  lightly  damped,  i.e.,  am  is  small,  it 

is  noted  that  in  the  neighborhood  of  w  *  a)m,  the  first 

term  in  equation  (3.25)  dominates,  and  Rm  is  nearly 

constant  for  small  changes  in  u.  Using  this  assump- 

til 

tion,  consider  the  nr  peak  in  the  magnitude  curve  of 
h(jw)  (see  Figure  3.11).  Pick  two  points  on  the  curve 
one  on  each  side  of  the  peak.  Now  then,  evaluate 


requency 


Figure  3.11.  Two  Point  Estimation 
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equation  (3.25)  at  each  of  the  two  points  to  get 


hl  - 

h(ju 

\  _  m 

i/ 

-  + 

sm 

(3. 

.26) 

h2  " 

h(joi 

i  -  ** 

2  Jw 2  ’ 

—  + 
sm 

Rm 

(3. 

.27) 

Solving  for 

a„  results  in 
m 

a  - 
m 

<hi 

-  Rm)(j“l  “ 

sm> 

- 

(h2 

-  V(j“2  “ 

sm> 

(3. 

.28) 

which  after  averaging  yields 

%  -  *  [<»1  -  V  U“1  - 

+  <h2  *  V<>2  '  sm>]  <3-29) 

Also,  sm  can  be  determined  from  equation  (3.28) 

>!<hl  -  V  '  >2(h2  '  Rm> 

sm - -"h2 -  <3'30) 

Since  these  equations  contain  complex  variables  we 
can  separate  them  into  real  and  imaginary  parts .  For 
convenience  in  evaluation,  define  the  following: 

a  -  ar  +  ja1 
m  m  J  m 

s  *  -a  +  1u 
m  m  J  m 


s 


h!  -  hl  +  jhj; 
h~  «  h~  +  jh~ 


(3.31 


Substituting  these  relationships  into  equations  (3.29) 
and  (3.30)  results  in: 

am  +  Jam  ’  *  [<hl  +  3hl  -  *£ 

-  J>£)<:Wl  +  V  -  J"m> 

+  <h|  +  jha  -  R* 

-  +  °m  -  J“m>]  <3-32 

.  *  Jh*  -  R*  -  JR*) 

-am  +  3aim  ”  hr  ,  ,.i  Tr  TTI 
h^  +  jh^  -  h2  _  jh£ 

_  j“2<h2  +  -  *£  - 

hl  +  3hl  “  h2  "  3h2  (3.3: 


Carrying  out  the  indicated  multiplications  and  divi¬ 
sions,  and  then  separating  into  real  and  imaginary 


parts  results  in: 


-(u>2  -  w1) 


m  (hT  -  ho)2  +  (ht  -  hh' 
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[Rm(hl  -  h2>  '  R>I  -  h2> 

+  h^h*  -  hjh![]  (3.34) 

„  _  1 

“m  (hj  -  +  (hj  - 

[<“2  -  “1>  *  h 
+  -  h^hj  + 

+  „2[<hj)2  -  h£h|  + 

«S  -  *[•.«>! +  *1  -  <> 

+  <.n<hJ  +  h|  -  2Ri) 

-  “2(hI  -  Rm>] 

am  ’  %[yhi  +  h|  -  2Ri) 

-  Vhl  +  h2  -  2Rm>  +  “l<hl  -  £> 

+“2(h2*Rm)]  <3-37> 

From  equations  (3.34)  through  (3.37)  it  can  be  seen 

that  sm  and  affl  can  be  found  if  Rm  is  known.  Recall 

that  R_  is  the  contribution  of  the  off-resonant  terms 
m 


I>  +  R,X  -  h§)] 

'  <^>2] 

<h2)2]]  (3.35) 

'  »i<ht  '  Rm> 

(3.36) 
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(equation  (3.24)).  If  an  estimate  is  made  for  R 

m 

then  values  for  s  and  a  could  be  obtained  and  then 

m  m 

an  iterative  technique  could  be  established  to  im¬ 
prove  the  estimates. 

The  technique  used  in  this  paper  is  to  as¬ 
sume  R  is  initially  zero  and  then  solve  for  s_  and 
am.  Then  using  these  values,  substitute  into  equa¬ 
tion  (3.24)  and  solve  for  a  new  Rm.  This  technique 
works  fine  when  all  of  the  modes  are  present  in  the 
frequency  range  of  analysis  and  modal  overlap  is 
light.  However,  in  real  structures  there  are  an  in¬ 
finite  number  of  modes,  all  of  which  cannot  be  mea¬ 
sured.  Therefore  it  may  be  necessary  to  consider  a 
contribution  to  the  transfer  functions  due  to  modes 
with  natural  frequencies  outside  of  the  frequency 
range  of  analysis.  This  contribution  can  be  con¬ 
sidered  either  to  be  a  complex  constant  or  a  complex 
function  of  frequency.  For  the  purpose  of  simplicity, 
we  will  consider  only  the  case  where  it  is  a  complex 
constant. 

The  constant  can  be  determined  by  a  least 
squares  fit  to  the  following: 
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h(jw) 


(3.38) 


where  Che  summation  includes  only  those  modes  in  the 
frequency  range  of  analysis.  Define  now  the  error 
function: 


e  -  h(jw) 


+ 


R 


(3.39) 


where  R  is  the  complex  least  squares  residual.  Then 
the  sum  of  the  error  squared  becomes: 


where 

k  *  block  size  of  the  frequency  response 
u)Q  ■  frequency  resolution  *  2irAf 


To  minimize  the  least  squares  residual  the  following 
must  be  true: 
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r  k/2  n 

r* i  .]  Li  ^  Z2kZ3k 

£{  (Z3k>"  +  <Z4k>: 


ZlkZ4k 


(Z,v)  +  <ZAU)‘ 


(k/2) 


(3.43) 


Now  then,  write  equation  (3.24)  to  include  the  con¬ 


stant  R  . 


x  XI 

[— 

-  sm  ]Z-ji  u<*>  - 


k  +  m 


*  T 

A. 

JO)  -  skJ 


(3.44) 


Again,  this  equation  can  be  separated  into  real  and 
imaginary  parts  for  convenience  of  computations .  De¬ 
fine  the  following: 


(<Jk)2  +  <U)k)2 
-2wo,_ 


(3.45) 
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Then  equation  (3.44)  can  be  written  as 


,  CTmam  +  am(“  +  Mm> 
(=0)2  +  <>»  +  »m)2 


'\_l  f  ZlkZ6k  ~  Z5kZ7kl 

Zj,  L  (Z6k^  +  <Z7k>"J 


+  R 


k  =  1 
k  ^  m 


m 


a  a*;  -  a*(u>  +  <u  ) 
mm  in _ m 

(Om)Z  +  (u  +  un)2 


V  r  Z5kZ6k  +  ZlkZ7k1  .  .'1 

.44  L<z6k)"  +  «7k>zJ 


k  #  m 


(3.46) 


The  curve  fitting  now  uses  equation  (3.46)  in 
lieu  of  equation  (3.24).  Again,  Rm  is  assumed  ini¬ 
tially  to  be  zero  and  then  values  for  s„  and  can 
be  computed.  Then  an  iterative  solution  can  be  ar¬ 
rived  at  using  equations  (3.34)  through  (3.37),  (3.43), 
and  (3.46).  A  computer  algorithm  using  these  equations 
can  be  found  in  Appendix  B. 

This  technique  works  fine  for  systems  with 
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light  to  moderate  damping.  When  there  is  heavy  modal 
overlap  due  to  closely  spaced  resonant  frequencies, 
but  with  light  damping,  the  technique  again  works  rea¬ 
sonably  well.  However  if  there  are  closely  spaced 
frequencies  and  the  damping  ratios  approach  0.2,  then 
the  error  in  the  estimated  parameters  begins  to  get 
too  large.  A  discussion  on  the  error  versus  damping 
ratio  is  presented  in  Chapter  IV  in  conjunction  with 
the  example  problems. 

3.3.2  Curve  Fitting  of  Quadrature  Response 

As  previously  mentioned,  the  total  response 
curve  fitting  technique  works  well  when  the  degree  of 
modal  overlap  is  light.  As  the  degree  of  overlap  in¬ 
creases  ,  this  method  yields  greater  errors .  One  meth¬ 
od  which  provides  more  accurate  results  is  to  curve 
fit  the  quadrature,  or  imaginary,  part  of  the  transfer 
function.  The  main  reason  for  choosing  the  quadrature 
response  is  the  rapid  change  it  displays  near  resonant 
frequencies.  The  quadrature  response  peaks  much 
sharper  at  a  resonance  as  compared  to  the  total  re¬ 
sponse  curve. 

In  this  section,  the  imaginary  curve  fit  algo¬ 
rithm  will  be  developed.  The  assumptions  used  appear 


58 


not  to  have  any  degrading  effects. 

Consider  the  summation  of  equation  (3.22). 

After  expanding  into  real  and  imaginary  parts  the  imag¬ 
inary  part  can  be  written  as : 


Due  to  the  form  of  the  above  equation,  a  closed  form 
solution  for  the  unknown  parameters  is  not  feasible. 
A  Newton  -  Raphson  technique  can  be  applied  in  the 
following  manner.  Define  the  following  function 
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estimates  were  available  for  these  parameters  then  an 
iterative  technique  could  be  applied  to  improve  the 
estimates.  It  can  be  seen  from  equation  (3.47)  that 


if  the  parameters  and  are  known  then  equation 
(3.47)  becomes  a  linear  function  of  and  8^. 
the  residues  (a^) .  Thus ,  it  would  only  be  necessary 
to  provide  initial  estimates  for  natural  frequencies 
and  damping  ratios.  The  estimates  for  the  residues 
(a^)  are  then  obtained  by  using  a  Gaussian  elimination 
technique  on  equation  (3.47)  evaluated  at  2n  frequency 
points  (preferably  two  points  on  each  peak) . 

In  order  to  define  an  iterative  technique  for 
the  solution  of  equation  (3.47),  the  function  g  (equa¬ 
tion  (3.48))  can  be  rewritten  as 


g  =  f(a 


__  j  | 

m  m  m 


V  +  f  (v 


U) 


k’ 


V  6k> 


where 


h1 

(3.49) 


f(a  ,  to  ,  a  ,  8  )  *  the  term  of  equation 
m  m  m  m 

(3.48)  associated  with  the 
mt^1  pole 

I 

f  (a^,  a^,  8^)  *  the  summation  of  all  terms 

in  equation  (3.48)  except 
for  the  term  associated 
with  the  m  pole 
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h^-  =  Im|h(jcL))J 

Using  initial  estimates  for  the  unknown  parameters, 
the  function  defined  by  equation  (3.49)  can  be  ex¬ 
panded  in  a  Taylor  series  expansion.  Neglecting  the 
second  and  higher  order  terms,  and  assuming  that  the 

*  th 

function  f  is  a  constant  at  frequencies  near  the  m 
resonant  frequency,  the  expansion  becomes: 


Sn+l  =  gn  +  Hi 
m 
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)  +  !s_ 

m _ .I  m  '  3u) 


(a  -  q 

nn+l  mn 
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(<*  -  a  >  +  T§~ 


m. 


n 


n+1  n 


m'  3  3. 
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(u> 

0) 

mn+l 

mn 

(B 

3  : 

mn+l 

“n 

(3.50) 


where  the  n  denotes  the  initial  estimate  and  n  +  1  is 
the  improved  estimate.  The  partial  derivatives  are 
evaluated  at  the  initial  estimates.  Since  gn  +  ^ 
should  equal  zero,  equation  (3.50)  can  be  written  as: 


Sn  3am 
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(Vi '  v> +  It 
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n  m 
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“n+l 


w  ) 
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(a  -  a  )  +  If-  (3  -  3  ) 

“n+l  mn  36m  I  ®n+l  mn 
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(3.51) 
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If  equation  (3.51)  is  written  at  four  frequency  points 
on  the  m  peak  then  a  matrix  equation  can  be  set  up 
as  follows: 


Ifi 

3gl 

3gl 

9a 

9  co 
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m 

m 

m 
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m 

m 
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Define  the  following  matrices : 


(3.52) 
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Using  these  definitions,  equation  (3.52)  can  be 
written  as: 


^n+1  ^n 


x,  ‘  J"18 

-n  —  2n 


(3.53) 


Thus,  with  initial  estimates  for  the  damping  ratios 
and  natural  frequencies,  equation  (3.47)  can  be  solved 
for  initial  estimates  of  the  residues.  Once  initial 
estimates  are  obtained  for  natural  frequencies,  damping 
ratios,  and  residues,  then  equation  (3.53)  can  be  used 
to  improve  the  estimates.  In  order  to  get  the  initial 
estimates  for  the  natural  frequencies  and  damping 
ratios  which  define  the  system  poles,  use  was  made  of 
the  method  of  Kennedy  and  Pancu  discussed  earlier. 

As  will  be  shown  later,  this  method  provides  fair 
estimates  for  the  system  poles,  but  not  necessarily 
the  mode  shapes.  For  all  cases  considered  in  this 
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thesis,  the  circle  fit  method  provided  estimates  with 
sufficient  accuracy  for  this  algorithm  to  converge. 

The  imaginary  curve  fitting  technique  provides 
very  accurate  estimates  for  the  modal  parameters  for 
systems  with  a  light  degree  of  modal  overlap.  As  the 
degree  of  modal  overlap  increases,  more  accurate  ini¬ 
tial  estimates  are  required  for  this  method  to  con¬ 
verge.  As  the  damping  ratio  approaches  0.2  with  reso¬ 
nant  frequencies  closely  spaced,  this  method  becomes 
unstable  if  the  initial  estimates  for  natural  fre¬ 
quency  and  damping  ratio  are  in  error  by  more  than 
107. . 

3.4  Summary 

Each  of  these  methods  has  its  advantages  and 
limitations.  Chapter  IV  will  present  two  example 
problems,  each  one  analyzed  by  each  curve  fit  method 
discussed  in  Sections  3.2.2,  3.3.1,  and  3.3.2.  ihe 
results  are  summarized  to  show  the  advantages  and 
limitations  of  each  method. 


CHAPTER  IV 


EXAMPLES  OF  CURVE  FITTING  TECHNIQUES 


4.1  Introduction 

This  chapter  will  be  devoted  to  exploring  the 
use  of  the  curve  fitting  techniques  described  in  Chap¬ 
ter  III  on  an  example  problem.  The  problem  to  be  con¬ 
sidered  will  be  the  two  degree  of  freedom  system  shown 
in  Figure  4.1. 


Figure  4.1.  Example  Problem 

This  system  will  be  used  for  two  example  problems,  one 
with  the  system  poles  separated  widely  and  one  with 
the  poles  close  together.  The  estimates  obtained  in 
each  case  will  be  compared  to  the  actual  results  in 


order  Co  show  the  usefulness  and  extent  of  each  curve 
fit  method. 


4.2  Example  Problem  One 

Define  the  variables  of  Figure  4.1  as: 


For  this  system,  the  actual  values  for  the  modal 
parameters  as  computed  via  a  complex  eigen-solver 
routine  are: 


U1 

C1 


0.97103  rad/sec 
0.00765 


“2 

?2 


10.29834  rad/ sec 
0.00980 


r 


mode  1 


0.05711  /-359.17 


1.0 


/  0C 


L 


1.0 


mode  2 


/  0C 


J 

1 


0.00964  / -171. 27c 


J 


Recall  that  the  transfer  matrix  can  be  written  as : 
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H(s) 


•u(s) 

hl2Cs) 

21(s) 

h22(s) 

If  the  frequency  response  functions  h^Cjw)  and  h2i(jw) 
are  measured,  then  all  of  the  modal  parameter  estima¬ 
tions  may  be  obtained.  In  all  of  the  cases  to  follow, 
h^  and  h2^  are  used  to  determine  the  estimates  of  the 
natural  frequencies,  damping  ratios,  and  mode  shapes. 


4.2.1  Method  of  Kennedy  and  Pancu  (Circle  Fit) 

The  algorithm  in  Appendix  A  was  used  to  compute 
the  estimations  shown  in  Table  4.1.  Immediately  the 


Table  4.1 


PARA¬ 

METER 


0>1 


Uif 


mode 

1 

mode 

2 


ESTIMATION 


0.96000 
0.03577 
10.29800 
0.00980 
0.05711  /-359 


/ 


.16°' 

o 


L 


0.00965  /-171.26 


o 


ACTUAL  VALUE 


0.97103 

0.00765 

10.29834 

0.00980 


0.05711  /-359.170' 
1  /  0° 


/ 


0.00964  /-171 


o 

.27° 


7.  ERROR 


1.14 

367.58 

0.003 

0.0 

0 . 0M*^ 
0.03P* 

0.11M 

0.01P 


M  *  Magnitude 


Phase 
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error  for  5^  stands  out.  This  gross  error  can  be 
readily  explained,  and  is  a  good  example  of  one  of  the 
problems  associated  with  this  method.  This  error  was 
the  result  of  insufficient  frequency  resolution  in  the 
frequency  response  data.  The  above  problem  was  re¬ 
worked  with  better  resolution  by  looking  only  at  the 
frequencies  near  u>^.  This  provided  a  0.01%  error  in 

Si- 

One  other  problem  associated  with  this  method 
is  that  of  one  mode  dominating  the  other  nearby  modes. 
For  this  problem,  even  though  mode  2  dominates,  it  is 
separated  far  enough  in  frequency  so  as  not  to  inter¬ 
act  with  mode  1.  Figure  4.2  is  the  polar  plot  for  the 
h^  frequency  response  of  this  problem.  It  appears  at 
first  to  be  a  plot  of  a  single  degree  of  freedom  sys¬ 
tem,  however,  a  closer  look  reveals  a  very  small  re¬ 
sponse  circle  near  the  origin.  This  domination  can 
be  seen  more  easily  on  the  log-magnitude  plot  of  Fig¬ 
ure  4.3.  This  type  of  domination  can  be  a  problem  if 
the  natural  frequencies  are  close  together. 

The  damping  ratio  in  this  problem  was  varied 
from  approximately  0.008  to  approximately  0.20.  The 
exact  solutions  are  tabulated  in  Appendix  D.l  and  the 
curve  fit  estimations  are  given  in  Appendix  D.2.  As 


a  summary  of  the  results,  the  maximum  errors  in  fre¬ 
quency,  damping,  and  mode  shape  estimations  are  plot¬ 
ted  in  Figure  4.4  as  a  function  of  damping  ratio.  In 


Damping  Ratio 


Figure  4.4.  Percent  Error  vs  Damping  Ratio 

Using  Circle  Fit,  Poles  Separated 
(w,  =  10.298  rad/ sec, 

0)2  *  0.971  rad/sec) 

general,  as  the  damping  increases,  so  do  the  errors 
in  the  estimations.  The  very  large  error  for  dis¬ 
cussed  earlier  (see  Table  4.1),  was  not  plotted  be¬ 
cause  of  the  misleading  effect  it  would  have  in  the 
error  plot. 

It  should  be  noted  that  the  error  is  dependent, 
not  only  on  the  damping,  but  also  on  the  frequency  res¬ 
olution,  as  was  seen  earlier.  As  the  frequency  reso- 


obtained  using  the  original  values  of  damping  as  de¬ 
fined  in  Section  4.2.  From  Table  4.2  it  can  be  seen 
that  this  method  provides  good  estimates  for  natural 
frequencies  and  damping  ratios.  However,  the  error 
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in  the  2nc*  node  shape  gets  larger  as  damping  increases . 
It  is  noted  that  the  algorithm  converged  very  slowly 
on  frequency  response  This  points  out  one  of  the 

problems  with  this  technique.  As  noted  earlier,  one 
mode  dominates  in  this  problem.  When  this  situation 
occurs,  the  initial  estimates  provided  by  the  two 
point  estimation  are  not  very  accurate.  As  a  result, 
the  convergence  is  either  slow  or  non-existent.  This 
same  problem  can  be  encountered  when  the  damping  ratio 
is  increased. 

As  in  the  previous  section,  the  damping  ratio 
was  varied  from  approximately  0.008  to  0.20  and  the 
results  tabulated  in  Appendix  D.3.  Figure  4.5  is  a 
plot  of  the  error  vs  damping  for  this  algorithm.  The 
error  in  the  mode  shapes  associated  with  the  lack  of 
convergence  were  not  plotted  in  order  to  avoid  mis¬ 
representation  of  the  error  vs  damping.  Again,  as 
the  damping  increased,  the  errors,  in  general,  in¬ 
creased  for  a  fixed  frequency  resolution. 

It  appears  that  the  main  problem  with  this 
method  is  that  the  initial  estimates  for  the  resi¬ 
dues  are  sometimes  not  good  enough  to  provide  con¬ 
vergence.  The  initial  estimates  for  the  residues  are 
calculated  using  the  original  assumption  that  off- 


resonant  inodes  do  not  contribute  to  the  transfer  func 
tion  at  any  given  resonance.  This  assumption  is  ques 
tionable  under  conditions  of  heavy  modal  overlap  or 


Figure  4.5.  Percent  Error  vs  Damping  Ratio 
Using  Complex  Curve  Fit,  Poles 
Separated  (cd^  *  10.298  rad/ sec, 
u>2  *  0.971  rad/sec) 

where  one  mode  is  dominant.  The  initial  residue  es¬ 
timates  could  possibly  be  improved  by  the  use  of  a 
technique  similar  to  that  used  in  the  next  section. 


4.2.3  Curve  Fit  of  Imaginary  Part  of  H(juQ 

The  algorithm  of  Appendix  C  was  used  to  com¬ 
pute  the  estimates  in  this  section.  The  convergence 
of  this  technique  is  highly  dependent  upon  high  grade 
initial  parameter  estimates.  For  this  algorithm,  es- 
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timates  are  required  for  the  natural  frequencies  and 
damping  ratios  only,  since  the  algorithm  computes  the 
initial  estimates  for  the  residues.  If  the  initial 
estimates  of  damping  and  natural  frequency  are  more 
than  10%  in  error,  the  algorithm  may  diverge  using 
this  technique.  For  the  purpose  of  this  thesis,  the 
initial  estimates  were  taken  from  the  method  of 
Kennedy  and  Pancu  (Circle  Fit,  Section  4.2.1).  Except 
for  the  case  of  inadequate  frequency  resolution,  the 
circle  fit  method  gave  good  results  for  natural  fre¬ 
quencies  and  damping  ratios.  The  results  obtained 
using  the  imaginary  curve  fit  to  improve  the  initial 
estimates  are  shown  in  Table  4.3. 


Table  4 . 3 


PARA¬ 

METER 

ESTIMATION 

ACTUAL  VALUE 

%  ERROR 

0)1 

0.97103 

0.97103 

0.0 

?1 

0.00765 

0.00765 

0.0 

0)2 

10.29334 

10.29834 

0.0 

0.00980 

0.00980 

0.0 

mode 

fo. 05712  / - 35 9 . 17° 

0.05711  /-359.170’ 

0.01M 

1 

o 

o 

H 

1  /  0° 

0.0P 

mode 

O 

O 

o 

o 

r*4 

0.0M 

2 

|o. 00964  /-171.290 

0.00964  /-171.270! 

0.01P 

_ 

L  J 

As  in  the  previous  sections,  the  damping  ratio 
was  varied  in  the  approximate  range  of  0.008  to  0.20 
and  the  results  tabulated  in  Appendix  D.4.  Figure  4.6 
plots  the  error  in  the  parameter  estimation  vs  damping 


U) 


Damping  Ratio 


Figure  4.6.  Percent  Error  vs  Damping  Ratio 

Using  Imaginary  Curve  Fit,  Poles 
Separated  (&^  *  10.298  rad/sec, 
a>2  =  0.971  rad/sec) 

ratio.  The  broken  line  shown  for  the  mode  shape  plot 
was  drawn  that  way  because  of  the  magnitudes  of  the 
error.  Notice  the  drastic  change  in  the  scale  for  7. 
error.  If  this  mode  shape  error  would  have  been  plot¬ 
ted  on  one  of  the  previous  plots  (Figures  4.4  or  4.5) 
it  would  have  appeared  as  a  straight  line  very  near 
0%  error.  In  fact,  if  the  error  criterion  were  based 
on  an  accuracy  of  only  three  decimal  places  instead  of 
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five  (as  used  in  this  study)  then  there  would  be  07o 
error  for  all  parameter  estimations  in  this  problem. 

The  damping  ratio  was  not  increased  any  further 
because  of  the  difficulty  in  obtaining  natural  fre¬ 
quencies  and  damping  ratios.  This  is  an  area  that 
needs  further  research  before  the  full  limits  of  the 
imaginary  curve  fit  algorithm  can  be  tested. 

It  can  be  seen  that  this  method  is  more  de¬ 
sirable  than  the  previous  two  methods,  especially  if 
the  mode  shapes  are  of  interest. 

4.3  Example  Problem  Two 

Using  the  same  system  as  described  in  Figure 
4.1,  define  the  following  variables: 

m^  *  3  ^2  *  1 

=  36  k2  ■  6 

c^  =  .15  C£  =  0.05 

For  this  system,  the  actual  values  for  the  modal  pa¬ 
rameters  as  computed  via  a  complex  eigen  solver  rou¬ 
tine  are: 

o)^  ■  2.16993  rad/sec  u)2  “  3.91039  rad/sec 
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C1  =  0.00764  52  =  0.01068 


mode  1 


mode  2  = 


0.64593 


/  0° 

/ -178. 94° 


~ 


The  main  difference  between  this  problem  and  problem 
one  is  spacing  of  the  system  poles.  For  this  problem 
the  poles  are  close  together  compared  to  the  separa¬ 
tion  in  problem  one. 

As  in  the  previous  example  problem,  the  fre¬ 
quency  response  functions  h-^  and  h2^  were  used  to  es¬ 
timate  the  modal  parameters .  The  frequency  resolution 
in  this  problem  was  approximately  0.01  rad/ sec. 


4.3.1  Method  of  Kennedy  and  Pancu  (Circle  Fit) 

As  in  Section  4.2.1,  the  estimates  for  the 
modal  parameters  were  computed  and  the  results  pre¬ 
sented  in  Table  4.4.  One  can  immediately  see  that  the 
estimates  are  very  good.  The  frequency  resolution  was 
sufficient  in  this  case  to  provide  good  estimates, 
with  the  damping  at  a  low  level.  However,  as  the 
amount  of  damping  increased,  the  accuracy  dropped  off 


rapidly.  The  damping  ratio  was  varied  over  the  ap¬ 
proximate  range  of  0.008  to  0.20  and  the  results 
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Table  4.4 


PARA¬ 

METER 


ESTIMATION 


ACTUAL  VALUE 


7o  ERROR 


U1 

51 

a)2 

^2 

mode 

1 

mode 

2 


2.17000 

0.00763 

3.91000 

0.01068 

[6.21510  /-359.41c 

.  1 


/ 


C 


/ 


0.64733  /-178 . 94 


2.16993 
0.00764 
3.91039 
0.01068 

0.21526  / -359 , 41c 
1  /  0° 


/ 


0.64593  /-178. 94c 


0.003 

0.13 

0.01 

0.0 

0.07M 

0.0P 

0.22M 

0.0P 


tabulated  in  Appendix  D.6.  For  purposes  of  comparison, 
the  exact  solutions  are  tabulated  in  Appendix  D.5. 
Figure  4.7  shows  how  the  errors  increased  with  damping. 
The  errors  in  this  problem  are  larger  than  those  in 
the  first  example  because  the  poles  are  closer  toget¬ 
her  and  the  modes  interact  more  with  each  other.  Fig¬ 
ure  4.8  shows  the  polar  plot  of  h^  with  a  damping 
ratio  of  approximately  0.008.  Figure  4.9  is  the  same 
frequency  response  function  except  that  the  damping 
ratio  is  now  approximately  0.20.  Figures  4.10  and 
4.11  are  the  log-magnitude  plots  of  the  functions  rep- 


79 


resented  in  Figures  4.8  and  4.9,  respectively.  The 


Figure  4.7.  Percent  Error  vs  Damping  Using 
Circle  Fit,  Poles  Close 
<u>,  »  2.170  rad/sec, 
u)2  *  3.910  rad/sec) 


interaction  of  the  two  modes  can  be  clearly  seen  in 
these  figures.  This  interaction  poses  a  problem  when 
using  this  technique,  because  the  responses  tend  to 
not  approximate  a  circle  at  points  removed  from  the 
natural  frequency.  This  causes  a  circle  of  improper 
dimensions  to  be  fit  through  the  data,  thus  giving 
rise  to  larger  errors  in  damping  ratios  and  mode 
shapes . 
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Figure  4.8.  Polar  Plot  of  h^(C  *  0.008) 


Magnitude  of  H(jui) 
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Figure  4.10.  Log-Magnitude  Plot 
of  huU  -  0.008) 


Figure  4.11.  Log-Magnitude  Plot  of 
huU  -  0.20) 


j 
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4.3.2  Complex  Curve  Fit 

The  estimates  using  this  technique  are  given  in 
Table  4.5  for  the  original  damping  defined  in  Section 
4.3.  As  the  damping  was  increased  from  approximately 


Table  4.5 


PARA¬ 

METER 

— 

ESTIMATION 

ACTUAL  VALUE 

Z  ERROR 

“l 

2.16993 

2.16993 

0.0 

*1 

0.00764 

0.00764 

0.0 

“2 

3.91038 

3.91039 

0.002 

^2 

0.01068 

0.01068 

0.0 

mode 

0.21533  / -359.41°"! 

10.21526  / -359.41°! 

1 

o 

o 

*“■4 

_ 

1  ' 

t-* 

o 

o 

mode 

1 

t— 1 

o 

o 

r  i  /  o° 

0.09M 

2 

0.64532  /-178 . 95°  j 

0.64593  /-178 . 94° 

0.01P 

0.008  to  approximately  0.20  the  errors  also  increased 
(see  Appendix  D.7).  The  plot  of  error  vs  damping  in 
Figure  4.12  shows  rapid  increase  in  error  as  the 
damping  increased.  As  discussed  in  Section  4.2.2,  the 
assumptions  used  in  this  method  are  questionable  with 
heavy  modal  overlap.  Figures  4.10  and  4.11  show  that 
there  is  considerable  overlap  between  the  two  modes, 
which  explains  the  large  increase  in  error  with  in¬ 
crease  in  the  damping  ratio. 


Damping  Ratio 


Figure  4.12.  Percent  Error  vs  Damping  Ratio 
Using  Complex  Curve  Fit,  Poles 
Close  (u>^  *2.17  rad/sec, 

<^2  *  3.91  rad/sec) 


4.3.3  Curve  Fit  of  Imaginary  Part  of  H(1to) 

The  use  of  the  algorithm  of  Appendix  C  pro¬ 
vided  the  results  shown  in  Table  4.6.  As  the  damping 
was  varied  from  0.008  to  0.20,  the  errors  did  not  in- 
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crease  as  in  the  previous  two  methods .  The  results 


Table  4.6 


PARA¬ 

METER 


mode 


mode 


ESTIMATION 

2.16994 

0.00754 

3.91038 

0.01031 


ACTUAL  VALUE 
2.16993 

0.00764 

3.91039 

0.01068 

.21526  / -359  A 
1  / _ 0° 


[0. 21535  /-359.370]  jo. 21526  /-359.410  1 

L  1  /  °°  j  L  1  /  °°  J 
[  1  /  0°  1  [  1  /  0°  1 
[o. 64594  /-178 . 87°j  [0. 64593  /-178.940 j 


%  ERROR 
0.0004 

1.29 

0.0003 

3.50 

0.04M 

0.01P 

0.002M 

0.04P 


for  the  various  amounts  of  damping  can  be  seen  in  Ap¬ 
pendix  D.8.  The  plot  of  error  vs  damping  in  Figure 
4.13  clearly  shows  the  superiority  of  this  method. 


Figure  4.13.  Percent  Error  vs  Damping  Ratio 

Using  Imaginary  Curve  Fit,  Poles 
Close  (u)-,-  2.17  rad/sec, 
u>2*  3 . 91irad/sec) 
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Caution  must  be  exercised  in  using  this  meth¬ 
od  and  choosing  the  initial  estimates  for  natural  fre¬ 
quency  and  damping  ratio.  If  the  initial  estimates  are 
greater  than  10%  in  error,  as  mentioned  earlier,  there 
is  a  possibility  that  the  method  will  diverge. 

4.4  Determination  of  Mass,  Stiffness, 
and  Damping  Matrices 

As  an  extension  of  example  problem  one,  the 
data  from  Section  4.2.3,  with  the  damping  ratio  approx¬ 
imately  equal  to  0.008,  was  used  to  determine  the  mass, 
stiffness,  and  damping  matrices.  This  example  will 
demonstrate  the  use  of  equations  (2.21),  (2.27),  and 
(2.31). 

It  has  been  shown  that  the  modal  properties 
can  be  determined  from  frequency  response  data  with 
relatively  good  accuracy  (depending  on  the  method 
used) .  The  residue  at  each  pole  is  measured  and  then 
used  to  determine  the  mode  shape.  Recall  from  equa¬ 
tion  (2.10)  that  the  residue  matrix  for  mode  k  can 
be  written  in  terms  of  the  normalized  mode  shape. 

4  *  ak^ 

Recall  that  the  residue  matrix  can  also  be  written  as 
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\  *  *k^k 

where  r^.  is  the  measured  residue  vector  of  mode  k 
normalized  in  accordance  with  equation  (2.17).  Thus 

* k&£  •  ^k 

If  we  again  normalize  the  residue  vector  in  the  fol¬ 
lowing  manner, 


we  have  in  essence  normalized  the  mode  shapes , 
i.e. 


-k  ,  _  ^k 


Thus  the  normalized  mode  shapes  needed  for  equations 
(2.21),  (2.27),  and  (2.31)  can  be  determined  from  the 
measured  residue  vectors.  This  normalization  was  per 
formed  on  the  measured  residues.  Since  the  residue 
matrix  is  measured  and  the  normalized  mode  shapes  ob¬ 
tained,  then  the  complex  scalers,  a^,  can  be  deter¬ 
mined  from  equation  (2.10).  Once  the  scalers  are 
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determined,  the  mass,  stiffness,  and  damping  matrices 
may  be  determined  using  equations  (2.21),  (2.27),  and 
(2.31). 


After  normalization,  the  data  from  Section 
4.2.3  was  used  to  predict  the  results  tabulated  in 
Table  4.7.  From  Table  4.7  it  can  be  seen  that  the 


Table  4.7 


MATRIX 


ESTIMATION 

ACTUAL  VALUE 

fl.  02  o.oo’ 

fi.oo  0.00 

[o.OO  6.01 

[0.00  6.00 

fl06 . 34  -6.05" 

[i.06.0  -6.0  ' 

[-6.05  6.01 

[-6.0  6.0_ 

To. 20  -0.10 

r 0.20  -0.10* 

[0.10  0.10 

[-0.10  0.10 

M 


estimations  for  the  matrices  are  very  good.  The  error 
in  these  estimations  will  be  dependent  on  how  well  the 
residue  matrix  can  be  measured. 

4.5  Summary 

This  completes  the  chapter  on  example  prob¬ 
lems.  These  problems  have  shown  the  limitations  of 
each  curve  fit  method,  and  also  have  verified  equa¬ 
tions  (2.21),  (2.27),  and  (2.31)  which  can  be  used 
to  determine  the  mass,  stiffness,  and  damping  ma- 
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trices.  The  next  chapter  will  review  the  results  pre¬ 
sented  and  draw  conclusions  on  which  curve  fit  method 
is  best  tinder  given  conditions. 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 


Four  techniques  were  described  in  Chapter  III 
for  estimating  the  natural  frequencies,  damping  ratios, 
and  mode  shapes  of  a  system  using  the  frequency  re¬ 
sponse  functions  measured  at  various  points  in  the  sys¬ 
tem.  The  first  method  discussed  dealt  only  with  the 
quadrature  response  and  is  only  useful  in  very  special 
cases.  The  other  three  techniques  were  investigated 
in  more  detail.  The  three  curve  fitting  methods  were 
used  to  solve  two  example  problems,  one  with  the  sys¬ 
tem  poles  separated  in  frequency  by  approximately  10 
rad/ sec,  and  the  other  with  the  poles  separated  in 
frequency  by  approximately  1.7  rad/ sec.  These  exam¬ 
ple  problems  pointed  out  the  usefulness  of  the  differ¬ 
ent  techniques. 

The  method  of  Kennedy  and  Pancu  (circle  fit) 
provided  very  good  estimates  for  natural  frequencies 
and  damping  ratios,  provided  that  the  frequency  reso¬ 
lution  is  adequate.  However,  if  the  mode  shapes  are 
also  desired,  this  method  begins  to  break  down.  It 
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was  shown  in  Sections  4.2.1  and  4.3.1  that,  when  there 
is  heavy  modal  overlap  or  when  one  mode  dominates,  the 
method  yields  large  errors .  For  example ,  with  the 
poles  separated  in  frequency  by  approximately  10  rad/ 
sec,  a  frequency  resolution  of  approximately  0.02  rad/ 
sec,  and  with  a  damping  ratio  of  approximately  0.008, 
the  error  in  the  damping  ratio  estimate  was  368%. 

This  was  due  to  inadequate  frequency  resolution  in  a 
problem  where  one  mode  dominated.  When  the  poles  were 
separated  by  approximately  1.7  rad/ sec  in  frequency 
and  with  the  damping  ratio  near  .20,  and  a  frequency 
resolution  of  approximately  0.01  rad/sec,  this  method 
produced  an  error  of  168%  in  the  mode  shape  of  mode  2. 
However,  under  the  same  conditions,  the  natural  fre¬ 
quency  estimations  were  in  error  by  less  than  3%  and 
the  damping  ratios  were  in  error  by  less  than  87.. 

It  is  concluded  that  the  method  of  Kennedy  and 
Pancu  is  reliable  for  estimating  the  natural  frequen¬ 
cies  and  damping  ratios  except  in  the  cases  where  one 
mode  dominates  or  when  the  damping  ratio  exceeds  0.20. 

The  complex  curve  fit  method  also  provided 
resonable  estimates  under  certain  conditions.  For 
instance,  when  the  poles  were  separated  in  frequency 
by  10  rad/sec,  with  8%  damping  and  0.02  rad/sec  fre- 
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quency  resolution,  the  errors  for  natural  frequencies, 
damping  ratios,  and  mode  shapes  were  all  less  than  17.. 

However,  with  the  damping  at  3.8%,  and  the  other  fac¬ 
tors  as  above,  this  method  produced  an  error  of  198% 
in  the  estimation  of  mode  shape  2  while  maintaining 
less  than  1%  error  in  other  parameters .  This  large 
error  was  due  to  a  lack  of  convergence  (in  50  itera¬ 
tions)  of  the  technique  used.  This  method  makes  use 
of  an  assumption  that  is  not  valid  when  one  mode  dom¬ 
inates  or  when  the  degree  of  modal  overlap  is  heavy. 

It  is  this  assumption  that  causes  the  large  errors  to 
be  generated.  For  example,  with  the  poles  separated 
in  frequency  by  1.7  rad/sec,  damping  set  at  0.20,  and 
a  frequency  resolution  of  0.01  rad/ sec,  this  method 
yielded  errors  in  mode  shapes  as  large  as  74%  and 
errors  in  natural  frequencies  and  damping  ratios  as 
high  as  16%  and  88%,  respectively. 

It  is  concluded  that  the  complex  curve  fit 
works  well  under  the  conditions  of  light  modal  over¬ 
lap  provided  that  no  one  mode  dominates  the  response. 

When  one  mode  dominates,  the  estimates  are  still  good 
for  the  natural  frequencies  and  damping  ratios,  but 
the  estimates  for  the  mode  shapes  are  not  reliable. 

When  the  poles  are  close  together  where  they  interact 
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to  form  the  total  response,  this  technique  is  not  very 
reliable.  The  errors,  again,  are  due  primarily  to 
poor  initial  estimates  (computed  based  on  an  assump¬ 
tion  that  is  invalid  for  heavy  modal  overlap)  for 
which  the  algorithm  converges  to  the  wrong  values. 

A  modified  complex  curve  fit  method  could  be 
investigated  by  using  a  method  similar  to  that  dis¬ 
cussed  in  Section  4.2.3  to  provide  the  initial  esti¬ 
mates  for  the  previously  described  complex  curve  fit. 
This  could  possibly  make  this  method  more  desireable. 

The  final  method  discussed,  that  of  curve  fit¬ 
ting  the  imaginary  part  of  the  frequency  response, 
provided  excellent  results  for  all  cases  considered. 
The  maximum  error  for  any  parameter  or  mode  shape  was 
less  than  4%  for  all  cases  considered.  These  cases 
included  the  poles  separated  in  frequency  by  10  rad / 
sec,  a  frequency  resolution  of  0.02  rad/sec,  and  the 
damping  ratio  varying  from  0.008  to  0.20.  They  also 
included  the  cases  where  the  poles  were  separated  in 
frequency  by  1.7  rad/sec  and  the  frequency  resolution 
was  0.01  rad/sec,  with  the  damping  ratio  varying  from 
0.008  to  0.20.  The  results  obtained  using  this  tech¬ 
nique  makes  this  method  the  most  desireable  of  the 
three  methods  discussed.  The  main  problem  associated 


with  this  method  is  that  it  diverges  if  the  initial  es 
timates  for  natural  frequencies  and  damping  ratios  are 
in  error  by  more  than  10%.  Therefore  caution  must  be 
used  when  determining  the  initial  estimates  required 
to  start  this  method.  It  has  been  found  that,  in  all 
cases  considered,  if  the  method  converged,  it  did  so 
with  excellent  results.  Therefore,  a  trial  and  error 
type  approach  could  be  taken,  if  necessary,  to  deter¬ 
mine  the  initial  estimates.  If  the  method  diverges, 
then  another  estimate  could  be  tried  until  convergence 
is  obtained. 

In  summary,  if  only  the  natural  frequencies 
and  damping  ratios  are  desired,  the  method  of  Kennedy 
and  Pancu  provides  acceptable  results,  except  where 
one  mode  is  dominant.  However,  if  the  mode  shapes  are 
also  desired,  then  the  curve  fit  of  the  imaginary  re¬ 
sponse  is  the  most  desireable  from  the  standpoint  of 
minimum  error. 
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APPENDIX  A 
CIRCLE  FIT  ALGORITHM 


SUBROUTINE  CIRCLE (XP , YP , WP , IPLT , NDOF) 

DIMENSION  WP(500) ,XP(500) ,YP(500) ,RX(5) ,RI(5) 
DIMENSION  S<3) ,WN(5) ,X(6) ,Y(6) ,W(6) ,Z(5) 

THIS  SUBROUTINE  SEARCHES  FOR  THE  MAXIMUM 
CHANGE  IN  ARC  LENGTH  ON  THE  POLAR 
PLOT  OF  THE  FREQUENCY  RESPONSE  FUNCTION. 

THE  FREQUENCY  POINT  AT  WHICH  THE  ARC 
LENGTH  IS  A  MAXIMUM  IS  DESIGNATED  AS 
THE  NATURAL  FREQUENCY.  A  CIRCLE  IS 
THEN  FIT  THROUGH  FIVE  POINTS;  THE 
NATURAL  FREQUENCY  POINT  AND  TWO  ADJ¬ 
ACENT  POINTS  ON  EACH  SIDE  OF  THE 
NATURAL  FREQUENCY.  TWO  OTHER  SUB¬ 
ROUTINES  ARE  USED  TO  DETERMINE  THE 
MAGNITUDE  AND  PHASE  OF  THE  MODE 
SHAPES  AND  THE  DAMPING  RATIOS. 

XP-THE  REAL  PART  OF  THE  FREQUENCY  RESPONSE 
YP-THE  IMAGINARY  PART  OF  THE  FREQUENCY  RESPONSE 
WP-THE  FREQUENCY  POINTS  AT  WHICH  THE  FREQUENCY 
RESPONSE  WAS  MEASURED 

IPLT-NUMBER  OF  POINTS  IN  THE  FREQUENCY  RESPONSE 
NDOF=THE  NUMBER  OF  DEGREES  OF  FREEDOM 
S-ARC  LENGTH 
WN=NATURAL  FREQUENCY 

X  AND  Y-THE  ARRAYS  OF  POINTS  WHICH  ARE  USED 
FOR  THE  CIRCLE  FIT 
Z-DAMPING  RATIO 

RX  AND  RI-THE  REAL  AND  IMAGINARY  PART  OF  THE 
MODE  SHAPE  COMPONENT,  RESPECTIVELY 

IC-0 

NUMPTS-IPLT-3 
DO  12  I-l.NUMPTS 
DO  1  J-1,3 

1  S (J) »SQRT( (ABS (XP (I+J) -XP(I+J- 1) ) )**2+ 

$  (ABS  (YP  ( I+J)  -YP  (I+J- 1)  ) )  **2) 

IF(S(l).LT.S(2).AND.S(3).LE.S(2))GO  TO  2 
GO  TO  10 

2  IC-IC+1 
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IF(IC.GT.NDOF)  PRINT  15 
IF(S(3) .GT.S(l))GO  TO  4 
WN(IC)-WP(I+1) 

DO  3  L-1,5 
W(L)»WP( I+L-2) 

X(L)=XP (I+L-2) 

3  Y(L)-YP (I+L-2) 

GO  TO  6 

4  WN(IC)«WP(I+2) 

DO  5  L-1,5 
W(L)*WP (I+L-l) 

X(L)»XP(I+L-1) 

5  Y(L)-YP( I+L-l) 

6  CONTINUE 

CALL  CIRCLE (X, Y.XCEN, YCEN, RADIUS) 

PRINT  7 

7  FORMAT (///,14X,*DATA  POINTS*, /, 5X,*OMEGA*, 

$4X , *REAL  PART  IMAG . PART* , / ) 

PRINT  8, (W(L) ,X(L) ,Y(L) , L-1,5) 

8  FORMAT ( 2X, F8 . 4 , 2F12 . 5) 

PRINT  9, XCEN, YCEN, RADIUS 

9  FORMAT ( / / / , 11X , *CIRCLE  PARAMETERS*, /, 5X, 

$*CENTER  (REAL , IMAG . ) * , 5X , *RADIUS* , / , 2X , 

$*(*  F10  5  *  *  FIO  5  *)*  F12  5) 

IF (ABS  ( Y ( 3 )  )  .  GT . ABS  (YCEN)  )  RADIUS-RADIUS 
CALL  DAMP (W ,  X ,  Y , XCEN , YCEN , RADIUS , ZETA , ALPHA) 
RX(IC)=RADIUS*COS (ALPHA) *2 . 

RI (IC) -RADIUS*SIN (ALPHA) *2 . 

Z(IC)-ZETA 

10  CONTINUE 

IF(I . EQ . NUMFTS . AND . IC . EQ. 0) PRINT  11 

11  FORMAT (///,5X,*N0  NATURAL  FREQUENCIES*,/, 
$5X,*NERE  FOUND*) 

12  CONTINUE 

IF (IC.LT.NDOF) PRINT  16 
IF (IC.EQ. 0)GO  TO  14 

PRINT  13, (I,WN(I) ,RX(I) ,RI(I) ,Z(I) , 1-1 , IC) 

13  FORMAT (/// , 5X,*WN(*, 12 , *)-*, FIO. 5 , 5X, 

$*MODE  SHAPE  COMPONENT  -(*,E13.5,*)  +  J(*. 
$E13 . 5 ,*)*,/, 5X, *DAMPING  RATIO-* , FIO . 5) 

14  CONTINUE 

15  FORMAT (///,5X,*THE  NUMBER  OF  NATURAL*,/, 

$*FREQUENCIES  FOUND  EXCEEDS  THE  NUMBER*, 

$/,*  OF  DEGREES  OF  FREEDOM*) 

16  FORMAT (///,5X,*THE  NUMBER  OF  NATURAL*,/, 

$*  FREQUENCIES  FOUND  IS  LESS  THAN  THE*,/, 

$*  NUMBER  OF  DEGREES  OF  FREEDOM*) 
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RETURN 

END 


SUBROUTINE  CIRFIT (X , Y , XCEN , YCEN , RADIUS) 
DIMENSION  X(5) , Y(5) ,A(3,3) ,F(3) ,R(3) 

REAL  ILL 

THIS  SUBROUTINE  CALCULATES  THE  CIRCLE  FIT 
PARAMETERS  WHICH  DETERMINE  THE  MAGNITUDE 
OF  THE  COMPLEX  RESIDUE. 

THE  VARIABLE  ILL  IS  USED  TO  DETERMINE  IF  THE 
MATRIX  IS  ILL-CONDITIONED. 

DO  1  1-1,3 
DO  1  J-1,3 
F(I)=0,0 
A(I,J)«0,0 
A(3 , 3)=5 , 0 
DO  2  1-1,5 

A(l, 1)-A(1, 1)+X(I)*X(I) 
A(1,2)-A(1,2)+X(I)*Y(I) 

A(1,3)-A(1,3)+X(I) 

A(2,2)»A(2,2)+Y(I)*Y(I) 

A(2,3)=A(2,3)+Y(I) 

F(1)*F(1) - (X(I)*X(I)+Y(I)*Y(I) )*X(I) 
F(2)»F(2)-(X(I)*X(I)+Y(I)*Y(I))*Y(I) 
F(3)-F(3)-(X(I)*X(I)+Y(I)*Y(I)) 

A(2 , 1)-A(1, 2) 

A(3,2)-A(2,3) 

A(3 , 1)-A(1, 3) 

ILL-(A(1,1)+2.*A(1,2)+2.*A(1,3)+A(2,2)  + 
$2.*A(2 , 3)+A(3 , 3)) /90000000 . 

DO  3  N-1,2 

DO  3  I-N+l , 3 

CONST-A(I,N)/A(N,N) 

IF (ABS (CONST) .LE.ABS( ILL))  PRINT  6 
F(I)«F(I)-F (N) *CONST 
DO  3  J-N+1,3 

A(I,J)-A(I,J)-A(N,J)*CONST 

CONTINUE 

R(3)-F(3)/A(3,3) 

DO  5  N-1,2 
M-3-N 

DO  4  I-M+l , 3 
F(M)-F(M)-A(M,I)*R(I) 

R(M)-F(M)/A(M,M) 

XCEN— R(l)/2. 


o  o  o  o  n  o  o 
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YCEN»-R(2)/2. 

RADIUS*SQRT (XCEN*XCEN+YCEN*YCEN- R(3 ) ) 

6  FORMAT (1H1,5X,*THE  COEFFICIENT  MATRIX  IS*,/, 

$5X, *ILL-CONDITIONED  OR  SINGULAR*, /, 5X, 

$*AND  THE  SOLUTION  MAY  BE  MEANINGLESS*) 

RETURN 

END 


SUBROUTINE  DAMP(W,X, Y,XC,YC,R,Z,A) 

DIMENSION  W(5) ,X(5) , Y(5) 

THIS  SUBROUTINE  DETERMINES  THE  DAMPING  RATIO 
AND  THE  PHASE  ANGLE  OF  THE  COMPLEX  RESIDUE. 
XC=X- COORDINATE  OF  THE  CIRCLE  CENTER 
YC*Y- COORDINATE  OF  THE  CIRCLE  CENTER 
R=RADIUS  OF  THE  CIRCLE 
Z=DAMPING  RATIO 
A=PHASE  ANGLE  IN  RADIANS 
CISQ=(X(3) -X(2) )*(X(3) -X(2))+(Y(3>- 
$Y(2)  )*(Y(3)  -Y(2)  ) 

C2SQ=  (X(4)  -  X(3)  )*(X(4)  -X(3)  )  +  (Y(4)  - 
$Y(3) )*(Y (4) -Y(3) ) 

Sl= (SQRT (ClSQ) ) / 2 . 

S2= (SQRT (C2SQ) ) / 2 . 

B1»SQRT(R**2-S1**2) 

B2=SQRT(R**2-S2**2) 

THETA1=2 . *ATAN2 (SI , Bl) 

THETA2-2 . *ATAN2(S2 , B2) 

PHI  1=*  (3 . 14159265-THETAI)  /2 . 

PHI2*(3 . 14159265+THETA2) /2 . 

BETA1*W(2)/W(3) 

BETA2=W(4)/W(3) 

ZI* ( ( 1 . -BETAI**2) *SIN (PHI1) / COS (PHI1) ) / 
$(2.*BETA1) 

Z2*( (1 . -BETA2**2)*SIN(PHI2) /COS (PHI2) ) / 
$(2.*BETA2) 

Z=(Zl+Z2)/2. 

XX2=ABS(X(3)-XC) 

YY2-ABS (Y(3)-YC) 

IF(Y(3) ,GT.YC)GO  TO  1 
A»ATAN2(XX2,YY2) 

IF(X(3)  .LT.XC)A—  A 
GO  TO  2 

1  A«ATAN2(YY2,XX2)+1. 57079632 
IF(X(3).LT.XC)A— A 

2  CONTINUE 


1 


1 
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RETURN 

END 
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APPENDIX  B 


MULTI  DEGREE  OF  FREEDOM  CURVE  FIT  ALGORITHM 


SUBROUTINE  HJW  (XP ,  YP ,  WP  , IPLT , NDOF) 

DIMENSION  WP(500) ,XP(500) ,YP(500) 

DIMENSION  W1 (5) ,W2(5) ,H1R(5) ,H1I(5) 
DIMENSION  H2R(5)  ,H2I(5)  , WN(5)  , RR(5)  ,RI(5) 
DIMENSION  SR(5 , 2) , SI(5) ,AR(5) ,AI(5) ,Z(5) 
THIS  SUBROUTINE  PICKS  OUT  THE  PEAKS  IN 
THE  IMAGINARY  PART  OF  THE  FREQUENCY 
RESPONSE.  TWO  POINTS,  ONE  ON  EACH 
SIDE  OF  THE  PEAK,  ARE  PICKED  FOR 
EACH  PEAK.  THESE  POINTS  ARE  THEN 
USED  IN  THE  TWO  POINT  ESTIMATION 
TECHNIQUE . 

XP-REAL  PART  OF  THE  FREQUENCY  RESPONSE 
YP- IMAGINARY  PART  OF  THE  FREQUENCY 
RESPONSE 

WP-FREQUENCY  POINTS  AT  WHICH  THE  FRE¬ 
QUENCY  RESPONSE  WAS  MEASURED 
IPLT=NUMBER  OF  POINTS  IN  THE  FREQUENCY 
RESPONSE 

NDOF=NUMBER  OF  DEGREES  OF  FREEDOM 
NUMPTS-IPT-2 
NPEAKS =0 
DO  6  I-l.NUMPTS 
IF (NPEAKS . GT . NDOF)  PRINT  29 
IF (ABS  (YP  (1+1)  )  .  GE .  ABS  (YP  (I)  )  .  AND .  ABS  (YP  (1+ 
$2)) .LE.ABS(YP(I+1)))  GO  TO  1 
GO  TO  4 

1  NPEAKS=NPEAKS+1 
DO  2  JJ-1,2 
IJK-I+l-JJ 
IJK-I+l+JJ 

RATIO 1-ABS  (YP (1+1)  /  YP  (IJK)  ) 

IF(RATIOl.GE. 1.0592537) GO  TO  3 

2  CONTINUE 

3  W1 (NPEAKS ) -WP ( IJK) 

W2 (NPEAKS )-WP( I JL) 

HIR (NPEAKS ) -XP ( I JK) 

H1I (NPEAKS) -YP (IJK) 
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H2R (NPEAKS ) =XP (IJL) 

H2I  (NPEAKS )  =YP  ( IJL) 

4  CONTINUE 

IF (NPEAKS. EQ.O. AND. I. EQ.NUMPTS) PRINT  5 

5  FORMAT (//,5X,*NO  NATURAL  FREQUENCIES* ,/, 5X , 
$*WERE  FOUND*) 

6  CONTINUE 

IF (NPEAKS. EQ.O) GO  TO  28 
PRINT  7, NPEAKS 

7  F0RMAT(1H1 , 5X, 12 , *  RESONANT  FREQUENCIES*,/, 
$5X, *WERE  FOUND*) 

PRINT  8 

8  FORMAT (//,5X,*POINTS  USED  IN  ESTIMATION*, 

$/ , 7X,*(2  PER  PEAK)*) 

DO  9  I J-l, NPEAKS 

9  PRINT  10 , IJ ,H1R(IJ) ,H1I (IJ) , IJ ,H2R(IJ) ,H2I(IJ) 

10  FORMAT (/,5X,*H1(*, 12,*)=  (*,E13.5,*)+ 

$J (*,E13 . 5 ,*)*,/ , 5X,*H2 (*, 12 ,*)=  (*, 
$E13.5,*)+J(*,E13.5,*)*) 

DO  11  I J-l, NPEAKS 
SR(IJ, 1)=0 . 0 
RR(IJ)=0 . 0 

11  RI(IJ)=0 . 0 
NN=1 

12  CONTINUE 

DO  13  MP=1, NPEAKS 
T=H1R(MP) -RR(MP) 

D=H1I(MP)-RI(MP) 

E=H2R(MP) -RR(MP) 

F-H2I(MP)-RI (MP) 

G=H1R(MP) -H2R(MP) 

P=H1I (MP) -H2I (MP) 

SR  (MP ,  2  )  =  (  (W2  (MP  )  *F  -  W1  (MP  )  *D)  *G+  (Wl  (MP  )  *T  - 
$W2 (MP) *E) *P) / (G*G+P*P) 

SI  (MP)  =  (  (Wl  (MP)  *T-W2  (MP) *E) *G-  (W2  (MP)*F- 
$W1 (MP) *D) *P ) / (G*G+P*P ) 

AR (MP ) =0 . 5* ( - ( SR (MP , 2 ) * (T+E)+D* (Wl (MP) - 
$SI(MP))+F*(W2 (MP) -SI (MP) ) ) ) 

AI (MP) =0 . 5*(T*(W1 (MP) -SI (MP) )+E*(W2 (MP) - 
$SI (MP) ) -SR(MP , 2) * (D+F) ) 

13  CONTINUE 
ICON-O 

DO  14  JCT-1, NPEAKS 

DIFF-ABS ( SR ( JCT , 2 ) - SR ( JCT , 1) ) 

IF(DIFF.LE. . 000001) ICON-ICON+1 
IF ( ICON. EQ. NPEAKS) GO  TO  25 

14  CONTINUE 
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RPI-0.0 

IF (NN.EQ. l)GO  TO  17 
DO  16  I-l.IPLT 
SUMR-0 . 0 
SUMI-0 . 0 

DO  15  J-l, NPEAKS 

Zl—  2.*(AR(J)*SR(J,2)+AI(J)*SI(J)) 
Z2«2.*WP(I)*AR(J) 

Z3-SR(J , 2)*SR(J, 2)+Sl (J)*SI(J) -WP(I)*WP(I) 
Z4-2.*WP(I)*SR(J,2) 

SUMR«SUMR+( (Z1*Z3-Z2*Z4)  /  (Z3*Z3+Z4*Z4)  ) 
SUMI-SUMI+( (Z2*Z3+Z1*Z4) / (Z3*Z3+Z4*Z4) ) 
CONTINUE 

RPR-RPR+XP ( I ) - SUMR 
RPI-RPI+YP (I) -SUMI 
CONTINUE 

RIPLT-FLOAT (IPLT) 

RPR=RPR/ RIPLT 
RP  I-RP  I  /  RIPLT 
CONTINUE 

DO  21  MR=1 , NPEAKS 
RBR=0 . 0 
RBI-0 . 0 
DO  20  NR-1,2 
WW—W1  (MR) 

IF (NR . EQ . 2 ) WW-W2 (MR) 

RBCR- (- (AR(MR) *SR(MR, 2) ) -AI (MR) *(WW+ 
$SI(MR) ) ) / (SR (MR, 2)*SR(MR, 2)+(WW+ 

$SI (MR) ) * (WW+SI (MR) ) ) 

RBCI- (AI (MR) *SR(MR,  2)-AR (MR) * (WW+ 

$SI  (MR)  )  )  /  (  SR  (MR ,  2  )  *SR  (MR ,  2  )  +  (WW+ 

$SI (MR) ) * (WW+SI (MR) ) ) 

SUMR-0 . 0 
SUMI-0 . 0 

IF (NPEAKS. EQ.l) TO  TO  19 

DO  18  KR-1, NPEAKS 

IF(KR.EQ.MR)GO  TO  18 

Ql—  2 .  *  (AR  (KR)  *SR  (KR ,  2  )  +AI  (KR)  *SI  (KR)  ) 

Q2-2 . *WW*AR(KR) 

Q3*SR(KR, 2)*SR(KR, 2)+SI(KR)*SI(KR) -WW*WW 
Q4-2 .*WW*SR(KR, 2) 

DENOM-Q  3  *Q3+Q4*Q4 
SUMR-SUMR+( (Q1*Q3-Q2*Q4) /DENOM) 

SUMI-SUMI+( (Q2*Q3+Q1*Q4) /DENOM) 

CONTINUE 

RBR-RBR+RBCR+SUMR 
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RBI=RBI+RBCI+SUMI 

20  CONTINUE 

RBR=RBR/2 . 

RBI-RBI/2. 

RR (MR) -RBR+RPR 

RI (MR) *RBI+RPI 

21  CONTINUE 

DO  22  IA=1 , NPEAKS 

22  SR(IA,1)«SR(IA,2) 

NN»NN+1 

IF (NN. GT . 50) TO  TO  23 

GO  TO  12 

23  PRINT  24.DIFF 

24  FORMAT (///  ,  5 X,  CONVERGENCE  NOT  OBTAINED*, 
$/ ,5X,*ON  THE  REAL  PART  OF  THE  POLES  ON*. 
$/,5X,*ALL  PEAKS  IN  50  ITERATIONS* ,//, 5X , 
$*DIFFERENCE  BETWEEN  LAST  TWO  ITERATIONS*, 
$/,5X,*ON  THE  REAL  PART  OF  THE  POLES*,/, 
$5X,*WAS  =*,F14. 6) 

25  CONTINUE 

IF (NPEAKS. LT.NDOF)  PRINT  30 

DO  26  LL»1, NPEAKS 

WN (LL) -SQRT (SR(LL , 2) *SR(LL , 2)+SI (LL) * 
$SI(LL) ) 

Z  (LL)  —  SR  (LL ,  2  )  /  WW  (LL) 

26  PRINT  27 ,LL,WN(LL) ,Z(LL) ,AR(LL) ,AI(LL) 

27  FORMAT (//,8X,*POLE  NUMBER  *,I2,/, 

$ 3X , *NATURAL  FREQUENCY  -*,F10.5,5X, 
$*DAMPING  RATIO  =*, F10 . 5 , / , 3X, 

$*MODE  SHAPE  COMPONENT  *  (*,E13.5, 

$*)  +/-  J(*,E13. 5,*)*) 

28  CONTINUE 

29  FORMAT (1H1,3X,*THE  NUMBER  OF  PEAKS*, /,3X, 
$*EXCEEDS  THE  NUMBER  OF  DEGREES  OF*,/,3X, 
$*FREEDOM*) 

30  FORMAT (//,3X,*THE  NUMBER  OF  PEAKS*, /,3X, 
$*IS  LESS  THAN  THE  NUMBER  OF  DEGREES*,/, 
$3X,*OF  FREEDOM*) 

RETURN 

END 
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APPENDIX  C 


IMAGINARY  CURVE  FIT  ALGORITHM 


SUBROUTINE  QUAD ( YP , WP , IPLT , NDOF , K , WN , Z) 
DIMENSION  WP(500) ,YP(500) ,F(10) ,U(10) 
DIMENSION  Wl(5) , W2 (5) , SNEW(5) ,W(10) 
DIMENSION  H3I(5) ,H4I(5) ,W3(5) ,W4(5) 
DIMENSION  H1I(5) ,Z(5) ,H2I(5) ,WN(5)  ’ 
DIMENSION  SR(5) ,SI(5) ,AR(5) ,AI(5) 

DIMENSION  ADD(4) ,AI(4) ,A2(4) 

REAL  JAC(5 ,10) 

THIS  SUBROUTINE  CURVE  FITS  THE  IMAG¬ 
INARY  PART  OF  THE  FREQUENCY  RESPONSE. 
INITIAL  ESTIMATES  FOR  THE  NATURAL 
FREQUENCIES  AND  DAMPING  RATIOS 
MUST  BE  INPUT. 

YP- IMAGINARY  PART  OF  THE  FREQUENCY 
RESPONSE 

WP-FREQUENCY  POINTS  AT  WHICH  THE 
RESPONSE  IS  MEASURED 
IPLT-NUMBER  OF  POINTS  IN  THE 
RESPONSE 

NDOF-NUMBER  OF  DEGREES  OF  FREEDOM 
K-NUMBER  OF  INITIAL  ESTIMATES  THAT 
ARE  INPUT 

WN-INITIAL  ESTIMATES  FOR  NATURAL 
FREQUENCIES 

Z-INITIAL  ESTIMATES  FOR  THE  DAMPING 
RATIOS 

NUMPTS-IPLT-2 
NPEAKS-0 
DO  4  I-l.NUMPTS 
IF (NPEAKS.GT. NDOF) PRINT  21 
IF(ABS(YP(I+1)) .GE.ABS(YP(I)) .AND.ABS(YP(I+ 
$2)) .LE.ABS(YP(I+l)))GO  TO  1 
GO  TO  2 
1  NPEAKS-NPEAKS+1 

W1 (NPEAKS ) =WP ( I - 7 ) 

W2(NPEAKS)-OT(I+9) 

HI I (NPEAKS) »YP (1-7) 

H2I (NPEAKS) -YP (1+9) 
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W3 (NPEAKS ) -WP ( I- 13 ) 

W4 (NPEAKS ) -WP ( 1+1 5 ) 

H3 I (NPEAKS ) - YP ( I - 1 3 ) 

H4 I (NPEAKS )*YP (1+15) 

WN (NPEAKS )-WP(I+l) 

2  CONTINUE 

IF (NPEAKS. EQ.O. AND. I. EQ.NUMPTS) PRINT  3 

3  FORMAT (//,5X,*NO  RESONANT  FREQUENCIES* , 

$1X, *WERE  FOUND*) 

4  CONTINUE 

IF (NPEAKS. EQ.O) GO  TO  22 
IF (NPEAKS. LT.NDOF)  PRINT  5 

5  FORMAT (1H1,5X,*THE  NUMBER  OF  PEAKS*, 

$1X,*F0UND  IS  LESS  THAN  THE  NUMBER*, 
$1X,*0R  DEGREES  OF  FREEDOM*) 

DO  6  I-l.K 
SR(I)— 2(I)*WN(I) 

6  SI(I)-WN(I)*SQRT(1.-Z(I)*Z(I)) 

DO  7  MP-1, NPEAKS 

IN-2*MP 

W(IN-l)-WKMP) 

W(IN)-W2(MP) 

F(IN-1)-H1I(MP) 

7  F(IN)-H2I(MP) 

CALL  GAUSS (SR, SI, NPEAKS, W,F,U) 

DO  8  MP-1, NPEAKS 
IN-2*MP 
SNEW(MP)-0 . 0 
AR(MP)-U(IN-1) 

8  AI(MP)-U(IN) 

PRINT  9 

9  FORMAT (1H1,3X,*INITIAL  VALUES*,//) 

DO  10  MP-1, NPEAKS 

10  PRINT  U,MP,WN(MP)  ,Z(MP)  ,AR(MP)  ,AI(MP) 

11  FORMAT (6X,*POLE  NO.*, 12, / , 2X, 

$*NATURAL  FREQUENCY-*, F10. 5, 5X, 

$*DAMPING  RATIO-* , F10 . 5 , / , 2X, 

$*RESIDUE  VECTOR-  (*,E13.5, 

$*)  */-  J(*,E13.5,*)*) 

NN-1 

12  CONTINUE 

DO  16  1-1, NPEAKS 
DO  15  NR-1,4 
ww-wi ( X  ^ 

IF(NR.EQ.2)NW»W2(I) 

IF (NR. EQ . 3)WW-W3 (I) 

IF (NR . EQ . 4) WW-W4 (I) 
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ADD(NR)*0. 

DO  14  KR-1, NPEAKS 
Ql= (2 . *AR(KR) *SI (KR) *SI (KR) - 
$2 . *AR (KR) *WW*WW- 2 . *SR (KR) *SR (KR) * 
$AR(KR) -4 . *AR(KR) *SI (KR)*AI (KR) ) *WW 
Q2- (SR(KR)*SR(KR)+2 . *SI (KR)*SI (KR)  + 

$2 . *WW*WW)* (SR(KR) *SR(KR) )+(SI (KR)* 

$SI (KR) ) * (SI (KR) *SI (KR) -2 . *WW*WW) +WW**4 
IF (NPEAKS . EQ . 1) GO  TO  13 
IF (KR.EQ. I) GO  TO  13 
ADD (NR) -ADD (NR) +Q1 / Q2 
GO  TO  14 

13  Zl(NR) «Q1 
Z2‘(NR)  -Q2 

14  CONTINUE 

JAC(NR,1)-(Z2(NR)*(-4.*SR(I)*AR(I)*WW- 
$4 . *SI (I) *AI (I) *WW) -Z1 (NR) * (4 . (SR(I)* 
$SR(I)*SR(I)+4.*SR(I)*SI(I)*SI(I)+4.* 
$SR(I)*WW*WW) ) / (Z2(NR) *Z2 (NR) ) 

JAC (NR, 2) - (Z2 (NR) * (4 . *WW*AR (I) *SI (I) - 
$4 . *WW*SR ( I ) *AI (I) ) -Z1 (NR) * (4 . *SR(I)* 
$SR(I)*SI (I)+4.*SI (I)*SI (I)*SI (I) - 
$4 . *SI (I) *WW*WW) ) / (Z2 (NR) *Z2 (NR) ) 
JAC(NR,3)=(2.*WW*(SI(I)*SI(I)-WW*WW- 
$SR(I)*SR(I)))/Z2(NR) 

JAC(NR,4)— 4.*SR(I)*SI(I)*WW/Z2(NR) 

15  CONTINUE 

G1=Z1(1)/Z2(1)4ADD(1)-H1I(I) 

G2=Z1(2) /Z2 (2)+ADD(2) -H2I(I) 
G3-Z1(3)/Z2(3)+ADD(3)-H3I(I) 

G4-Z1(4) / Z2 (4)+ADD(4) -H4I(I) 

ISIZE-4 

CALL  INVERT ( JAC , DMT , ISIZE ,15,5) 
SR(I)»SR(I)«(JAC(1,1)*G1+JAC(1,2)*G2* 
$JAC(1,3)*G3+JAC(1,4)*G4) 

SI (I) «SI (I) - (JAC (2 , 1)*G1+JAC(2 , 2)*G2+ 

$ JAC (2,3) *G3+JAC ( 2 , 4) *G4) 

AR(I)-AR(I) - (JAC (3 , 1)*G1+JAC(3 , 2)*G2+ 
$JAC(3 , 3)*G3+JAC(3 , 4)*G4) 

AI  (I) -AI  (I)  -  (JAC  (4 , 1) *G1+JAC  (4 , 2) *G2+ 

$ JAC (4,3) *G3+JAC (4 , 4) *G4) 

16  CONTINUE 
ICON-O 

DO  17  J-l, NPEAKS 
DIFF-ABS (SR(J) -SNEW(J) ) 

IF (DIFF . LE . . 000001) ICON-ICON+1 
IF (ICON. EQ. NPEAKS)  GO  TO  18 
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17  SNEW(J)=SR(J) 

NN-NN+1 

IF (NN. GT . 50) GO  TO  18 
GO  TO  12 

18  CONTINUE 
PRINT  19 

19  FORMAT (/ till,  3X, *FINAL  VALUES*,//) 

DO  20  I-l.NPEAKS 

WN(I)-SQRT(SR(I)*SR(I)+SI(I)*SI(I)) 
Z(I)  —  SR(I)/WN(I) 

20  PRINT  11,I,WN(I) ,Z(I) , AR(I) , AI (I) 

21  FORMAT (1H1,3X,*THE  NUMBER  OF  PEAKS*, 

$1X,*F0UND  EXCEEDS  THE  NUMBER  OF*, 
$1X,*DEGREES  OF  FREEDOM*) 

22  CONTINUE 
RETURN 
END 


SUBROUTINE  INVERT (A , DET , N , NDIGIT , NDIM) 
DIMENSION  A (NDIM, 8) 

C  THIS  SUBROUTINE  INVERTS  A  MATRIX 

C  AND  CALCULATES  THE  DETERMINANT. 

SUM-0. 

DO  1  1=1, N 
DO  1  J-l.N 

1  SUM-SUM+ABS (A( I , J) ) 

SUM-10 . ** (-NDIGIT/ 2 . )*SUM/N**2 

NP1-N+1 

NPN-N+N 

DO  2  1=1 ,N 

IPN-l+N 

DO  2  J-NP1 , NPN 
A(I, J)*0 . 

2  IF(IPN.EQ.J)  A(I, J)-l. 

DET-1. 

INTCH-0 
DO  9  I-1,N 
IP1-I+1 

IF(I.EQ.N)  GO  TO  5 
M-I 

DO  3  J-IPl.N 

3  IF(ABS(A(M,I)) .LT . ABS (A(J , I) ) )M-J 
IF(M.EQ.I)  GO  TO  5 
INTCH-INTCH+1 

DO  4  J-l.NPN 
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TEMP=A(M,J) 

A(M,J)-A(I,J) 

4  A(I , J) =TEMP 

5  IF (A(I ,1) . EQ. 0 . )  GO  TO  11 

IF (ABS(A(I , I) ) . LT . SUM)  PRINT  15 
DO  6  J-IPl.NPN 

6  A(I,J)=A(I,J)/A(I,I) 

DO  8  J-lfN 
IF(J.EQ.l)  GO  TO  8 
DO  7  K-IPl.NPN 

7  A(J,K)=A(J,K)-A(J,I)*A(I,K) 

A(J,I)=0. 

8  CONTINUE 

9  DET=DET*A(I ,  I) 

DET= (-1) **INTCH*DET 
DO  10  1=1, N 

DO  10  J=1,N 

10  A(I, J)=A(I, J+N) 

GO  TO  13 

11  PRINT  14 

DO  12  1=1, N 
DO  12  J=1,N 

12  A(I,J)»1. 

13  CONTINUE 

14  FORMAT (//,10X,*THE  MATRIX  IS  SINGULAR*, 
$1X,*AND  THEREFORE  NO  INVERSE  EXISTS*, 

$/ , 10X,*THE  ELEMENTS  OF  THE  MATRIX*, IX, 
$*ARE  ARBITRARILY  SET  TO  1*,/) 

15  FORMAT (//,10X,*THE  MATRIX  IS  ILL-CONDI*, 
$*TIONED  OR  SINGULAR  AND  THE  SOLUTION*. 
$1X, *MAY  BE  MEANINGLESS*,/) 

RETURN 

END 


SUBROUTINE  GAUS S  (X ,  Y ,  M ,  W ,  F ,  U) 

DIMENSION  X(l) , Y(l) ,W(1) ,F(1) ,U(1) ,D(10, 10) 

N2=2*M 

DO  1  N=1,M 

I=2*N 

DO  1  K=1,N2 

TOP 1=2 , *W(K) * (Y (N) *Y (N) -X(N) *X (N) - 
$W(K)*W(K)  ) 

TOP  2— 4 .  *W  (K)  *X  (N)  *Y  (N) 

BOT=  (X  (N)  *X  (N)  )  **2+(Y  (N)  *Y  (N)  )  **2+ 

$  (W(K)  *W(K)  )  **2+2 .  *X  (N)  *X  (N)  *Y  (N)  * 
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$Y(N)+2.*W(K)*W(K)*(X(N)*X(N)- 

$Y(N)*Y(N)) 

D(K,I-1)-TOP1/BOT 

D(K,I)»TOP2/BOT 

CONTINUE 

N3-N2-1 

DO  3  N-1.N3 

DO  3  I-N+1.N3 

CONST-D(I,N)/D(N,N) 

IF (ABS (CONST) .LE. . 0000000001)PRINT  2 
FORMAT (//,3X,*THE  D-MATRIX  IS  ILL-CONDI* 
$*TIONED*) 

F(I)»F(I)-F(N)*CONST 
DO  3  J-N+1.N2 

D(I,J)*D(I,J) -D(N, J)*CONST 
CONTINUE 

U(N2)-F(N2)/D(N2,N2) 

DO  5  NI-1.N3 

N-N2-NI 

DO  4  I-N+1.N2 

F(N)-F(N)-D(N,I)*U(I) 

U(N)-F(N)/D(N,N) 

RETURN 

END 


APPENDIX  D 


DATA 


D. 1  Exact  Solutions  -  Poles  Separated 


MULTIPLIER  OF  DAMPING  MATRIX 


PARAMETER 


ACTUAL  VALUE 


“1 

Cl 

^2 

^2 

mode 

1 

mode 

2 


0.97103 

0.00765 

10.29834 

0.00980 

'0.05711  /-359.17 

.  i  /  p— 
1  /  o 

0.00964  FT7T.T7 


r. 


MULTIPLIER  OF  DAMPING  MATRIX  -  5  ? 

PARAMETER 

ACTUAL  VALUE 

Ui 


0)/ 


0.97118 

0.03827 

10.29677 


C 


2 


0.04900 


mode 

1 

mode 

2 


0.05714  /-3S5.85 

1  j  0 == 

1  /  0 

0.01225  7-143.37 

k  ■ 
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I  MULTIPLIER  OF  DAMPING  MATRIX  -  16  1 

PARAMETER 

ACTUAL  VALUE 

*1 

0.97164 

0.07657 

“2 

10.29186 

^2 

0.09803 

mode 

TO. 05711  / -351 . 70 

1 

1  7  0 

mode 

1/0 

2 

0.01812  7^125756 

1  MULTIPLIER  OF  DAMPING  MATRIX  *  15 - 1 

PARAMETER 

ACTUAL  VALUE 

W1 

0.97242 

0.11493 

"2 

10.28362 

h 

0.14715 

mode 

TO. 05735  /-347 . 52  “ 

1 

1  1  7 - 0 - 

mode 

r  1  /  0  i 

2 

[0.02507  7-119 .'42 

(”  MULTIPLIER  OF  DAMPlNC  MATRIX  *  20 - 1 

PARAMETER 

ACTUAL  VALUE 

"1 

0.97352 

C! 

0.15338 

“2 

10.27200 

^2 

0.19639 

mode 

’0.05755  /-343.30  1 

1 

1  7  '  0 

mode 

1/01 

2 

0.03245  7-116.59 

1 - 

D.2  Circle  Fit  Data  -  Poles  Separated 


MULTIPLIER  OF  DAMPING  MATRIX 


PARAMETER 

ESTIMATION 

%  ERROR 

“i 

0.96000 

1.14 

0.03577 

367.58 

Cl)  ^ 

10.29800 

0.003 

c2 

0.00980 

0.0 

mode 

| 

TO. 05711  /-359.161 

0.0M* 

1 

i  1  7 - 0 - 

0.003P 

mode 

r  i/oi 

0.11M 

2 

[0.00965  7-T7T.26 

0.01P 

*M  *  Magnitude,  P 


PARAMETER 


mode 

1 

mode 

2 


ING  MATR 


ESTIMATION 


0.96000 

0.04185 

10.27500 

0.04915 


"0.05704  / - 355 . 8l" 

i  /  d 

:  i  /  o  : 

0.01247  7-iyy.ig 


1  ERROR 


1.15 

9.35 

0.21 

0.31 

0.17M 

0.01P 

1.81M 

0.09P 


D.3  Complex  Curve  Fit  -  Poles  Separated 


PARAMETER 


Wi 


mode 

1 

mode 

2 


ESTIMATION 


0.97201 

0.07664 

10.29230 

0.09778 

0.05723  / -352 . 06 

1  /  ~P 

1 

0.01822 


%  ERROR 


0.04 

0.09 

0.004 

0.26 

0.02M 

0.10P 

0.56M 

0.79P 


aaBiiaBiaagaEagaEaBaiaE 


PARAMETER 


“1 

*1 

(Do 


mode 

1 


mode 

2 


PARAMETER 


"1 

h 

(Do 


mode 

1 


mode 

2 


ESTIMATION 


0.98066 

0.15382 

10.28007 

0.19437 

0.05731  /-346.91 

1  Z  B  I 

1 

0.03573 


%  ERROR 


0.73 

0.29 

0.08 

1.03 

0.42M 

1.05P 

10 . 11M 
1.63P 


D.4  Imaginary  Curve  Fit  -  Poles  Separated 


PARAMETER 


mode 

1 

mode 

2 


ESTIMATION 

0.97164 

0.07657 

10.29186 

0.09803 

,05722  /-351. 

1  7  C 


X  ERROR 


722  /-351.70' 
1  0  ' 

1  /  0  = 

.01812  /-I26.4  ' 


0.0M 

0.0P 

0.0M 

0.0P 


PARAMETER 


mode 

1 

mode 

2 


ESTIMATION 

0.97242 

0.11493 

10.28362 

0.14715 

.05735  /- 347.52' 

1  7 - 0 - 


7.  ERROR 


1  / 
02507  T 


MULTIPLIER 
PARAMETER  I 


)F  DAMPING  MAT! 
ESTIMATION 
0.97352 


mode 

1 

mode 

2 


0.15338 
10.27200 
0.19639 
,05754  /- 343, 

i  Pnr 


03245  TZ 


0.02M 

0.0P 

0.0M 

0.0P 


X  ERROR 

0.0 


0.02M 

0.0P 

0.0M 

0.01P 


Circle  Fit  Data  -  Poles  Close 


PARAMETER 


Wi 


mode 

1 

mode 

2 


ESTIMATION 


2.16000 

0.07550 

3.89000 

0.10655 

0.20627  /-345 .76 
1 


1 

0.81231 


rm*im 


7.  ERROR 


0.58 

1.07 

0.40 

0.34 

4.82M 

2.36P 

22.42M 

3.70P 


PARAMETER 


(U-i 


mode 

1 

mode 

2 


ESTIMATION 


2.15000 

0.11490 

3.87000 

0.16058 

”0.20783  /- 342.45" 

1  T  TF 

1  /  0 

1.07109  7-146.85 


7.  ERROR 


1.20 

0.43 

0.75 

0.02 

4.92M 

2.49P 

56.18M 

10.83P 


PLIER 


D.7  Complex  Curve  Fit  -  Poles  Close 


tMufciaMaaajg 


PARAMETER 


U-i 


mode 

1 

mode 

2 


ESTIMATION  %  ERROR 


2.16993  0.0 

0.00764  0.0 

3.91038  0.0002 

0.01068  0.0 

0.21533  /-359.411  0.03M 

1  /  0  J  0.0P 

1  /  0  “|  0.09M 

0.64532  7-178.95  0.01P 


R  OF  DAMPING  MAT 


PARAMETER 


Wi 


PARAMETER 


(Di 


TRIX 


sisaaiaaiaaEaKasa 


ING  MATRIX 


ESTIMATION 


2.19491 

0.13419 

3.83325 

0.17129 


0.27998  / -356. 

1  7 - P" 

1  / _ 0 

0.52214  f- I6T 


.  02  j 


7.  ERROR 


0.86 

17.29 

1.69 

6.69 

28 . 09M 
1.38P 

23 . 86M 
0.36P 


ESTIMATION 


2.52032 

0.12580 

4.00429 

0.40395 

0.27337  /  269.06 

1  / 

1 

0.20754 


7.  ERROR 


15.54 

17.46 

2.94 

88.36 

23 . 54M 
73.89P 

71.06M 

41.16P 


MULTIPLIER 


PARAMETER 


OF  DAMPING  MATRIX 


ESTIMATION 


2.17306 

0.07538 

3.90476 

0.10317 


mode 

1 

mode 

2 


21677  / - 353 

1  7 - g 

1 

,66462 


=] 


%  ERROR 


0.02 

1.23 

0.02 

3.50 

i.03M 
M1P 

M6M 
0.41P 


Eiijmwreai 


PARAMETER 


OF  DAMPING  MAT 


ESTIMATION 


2.17708 


mode 

1 

mode 

2 


0.11301 

3.89756 

0.15496 

21872  / -350 . 59" 
1  /  "0~ 

1  /  o 

,68820  / -163  .'69 


7o  ERROR 


0.04 

1.23 

0.04 

3.48 


0.61P 


PARAMETER 


mode 

1 

mode 

2 


ESTIMATION 


2.18285 

0.15055 

3.88725 

0.20701 


,22154  /- 

i  r 


/-347.451 

z=rzj 


1  / _ 0 

72128  /-158 


%  ERROR 


0.07 

1.22 

0.07 

3.47 


0.82P 
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